STATISTICAL  ERROR  ANALYSIS  OF  A  DOA  ESTIMATOR  FOR  A 
PCL  SYSTEM  USING  THE  CRAMER-RAO  BOUND  THEOREM 

THESIS 


Kerim  SAY,  First  Lieutenant,  TUAF 


AFIT/GE/ENG/02M-23 


DEPARTMENT  OF  THE  AIR  FORCE 
AIR  UNIVERSITY 

AIR  FORCE  INSTITUTE  OF  TECHNOLOGY 


Wright-Patterson  Air  Force  Base,  Ohio 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 


Report  Classification 

Classification  of  this  page 

unclassified 

unclassified 

Classification  of  Abstract 

Limitation  of  Abstract 

unclassified 

UU 

Number  of  Pages 

111 

The  views  expressed  in  this  document  are  those  of  the  author  and  do  not 
reflect  the  official  policy  or  position  of  the  United  States  Air  Force,  Department  of 
Defense,  United  States  Government,  the  corresponding  agencies  of  any  other 
government,  NATO,  or  any  other  defense  organization. 

This  document  represents  the  results  of  research  based  on  information 
obtained  solely  from  open  sources.  No  agency,  whether  United  States 
Government  or  otherwise,  provided  any  threat  system  parameters,  or  weapons 
systems  performance  data  in  support  of  the  research  documented  herein. 


AFIT/GE/ENG/02M-23 


STATISTICAL  ERROR  ANALYSIS  OF  A  DOA  ESTIMATOR  FOR  A  PCL 
SYSTEM  USING  THE  CRAMER-RAO  BOUND  THEOREM 


THESIS 


Presented  to  the  Faculty 

Graduate  School  of  Engineering  and  Management 
Air  Force  Institute  of  Technology 
Air  University 

Air  Education  and  Training  Command 
In  Partial  Fulfillment  of  the  Requirements  for  the 
Degree  of  Master  of  Science  in  Electrical  Engineering 

Kerim  SAY,  B.S. 

First  Lieutenant,  TUAF 

March  2002 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 


AFIT/GE/ENG/02M-23 


STATISTICAL  ERROR  ANALYSIS  OF  A  DOA  ESTIMATOR  FOR  A 
PCL  SYSTEM  USING  THE  CRAMER-RAO  BOUND  THEOREM 


Kerim  SAY,  B.S. 
First  Lieutenant,  TUAF 


Approved: 

Maj.  Roger  L.  Claypoole,  Jr.,  PhD. 
Committee  Chairman 


IS  MAje  02_ 
Date 


Date 


16 

Date 


iii 


Acknowledgements 

I  would  like  to  thank  my  thesis  advisor  Maj.  Roger  L.  Claypoole,  Jr.  for 
his  support  and  guidance  throughout  this  research.  I  am  grateful  for  the 
friendly  working  environment  that  he  created.  The  contributions  of  my 
committee  Dr.  Andrew  J.  Terzuoli  and  Dr.  Michael  A.  Temple  were  also 
invaluable  and  certainly  appreciated.  I  especially  wish  to  thank  my  sponsor 
Dr.  Paul  E.  Howland,  who  made  contributions  in  the  PCL  area,  for  his 
support  from  great  distances.  I  also  would  like  to  thank  Dr.  Erdogan 
Dilaveroglu  from  Uludag  University,  Bursa/ Turkey  and  Mr.  Nicholas  J. 
Willis  for  their  advice  and  support. 

The  support  and  love  of  my  family  and  my  girlfriend  were  one  of  the 
major  keystones  in  completing  this  research.  Despite  the  great  geographical 
distances,  they  were  always  there  for  me.  I  have  enjoyed  the  sincere  and 
professional  working  environment  in  AEIT.  Therefore,  I  would  like  to  extend 
my  thanks  to  my  colleagues  from  the  U.S.  Military  Services,  to  the  family  of 
international  officers,  and  especially  to  the  Turkish  Air  Eorce  officers  and 
their  families  for  keeping  the  spirit  and  motivation  high.  I  want  to  take  the 
opportunity  to  thank  Lt.  Ahmet  Ozgetin  and  Lt.  Bari§  ^ahkoglu  for  their 
professionalism  and  true  friendship  during  our  study  sessions.  I  cannot 
finish  without  thanking  all  my  true  friends  both  here  and  back  in  Turkey; 
their  friendship  is  and  always  will  be  invaluable  to  me. 


IV 


Table  of  Contents 


Acknowledgements . iv 

Table  of  Contents . v 

List  of  Figures . ix 

Notation . xii 

Acronym  List . xiii 

Abstract . xiv 


STATISTICAL  ERROR  ANALYSIS  OF  A  DOA  ESTIMATOR  FOR  A  PCL 


SYSTEM  USING  THE  CRAMER-RAO  BOUND  THEOREM . 1 

1 .  INTRODUCTION . 1 

1.1  Background . 1 

1.2  Problem  Statement . 2 

1.3  Scope . 7 

1.4  Approach/ Methodology . 8 

1.5  System  Description . 9 

1.6  Assumptions . 10 

2  .  BACKGROUND . 12 

2.1  Development  of  RADAR  and  Basic  Terms . 12 


V 


2.2  Passive  Coherent  Location  System . 17 

2.2.1  Types  of  PCL . 18 

2.2.2  System  Operation . 18 

2.2.3  Characteristics  of  PCL . 19 

2.2.4  Advantages  of  PCL . 21 

2.2.5  Disadvantages  of  PCL . 22 

2.3  Cramer-Rao  Lower  Bound  (CRLB,  CRB) . 24 

2.3.1  Unbiased  Estimators . 25 

2.3.2  Minimum  Variance  Unbiased  (MVUB)  Estimator . 25 

2.3.3  Likelihood,  Log-likelihood,  and  Score  Eunctions . 27 

2.3.4  Cramer-Rao  Bound  Theorem . 29 

2.4  Phased  Array  Antennas . 31 

2.5  Summary . 33 

3  .  METHODOLOGY . 34 

3.1  System  Model . 34 

3.2  Derivation  of  CRB  for  Arbitrary  Signals  (CRB as) . 35 

3.3  Derivation  of  the  CRB  for  Constant  Modulus  Signals  (CRBcm) . 45 

3.4  CRB  Calculation . 50 

3.5  Summary . 52 

4  .  CONTRIBUTIONS  OF  RESEARCH . 53 

4.1  Scenario  1:  Single  Signal . 53 

4.2  Scenario  2:  Two  Signal  Sources . 55 

vi 


4.2.1  CRB  AS  of  Two  Signal  Sources . 56 

4.2.2  CRBcm  of  Two  Signal  Sources . 60 

4.3  Scenario  3:  Two  Signal  Sources  and  a  Transmitter  of  Opportunity 

(Interferer) . 67 

4.3.1  CRB  AS  of  Two  Signal  Sources  and  a  Transmitter  of  Opportunity 

(Interferer) . 67 

4.3.2  CRBcm  of  Two  Signal  Sources  and  a  Transmitter  of  Opportunity 

(Interferer) . 69 

4.4  Scenario  4:  Moving  Object  vs.  Different  Transmitter  of  Opportunity 

(Interferer)  Locations . 72 

4.4.1  CRB  AS  of  the  Moving  Object  for  Different  Transmitter  Locations72 

4.4.2  CRBcm  of  the  Moving  Object  for  Different  Transmitter  Locations76 

5  .  DISCUSSION . 80 

5.1  Conclusions . 80 

5.2  Future  Work . 81 

Appendix  A . 84 

Appendix  B . 86 

Appendix  C . 87 

Appendix  D . 88 


Vll 


Appendix  E 


90 


Vlll 


List  of  Figures 


Figure  1.  Transmit  and  Receive  Angles  (0t,  0r) . 3 

Figure  2.  Isorange  Ellipsoid . 3 


Figure  3.  Isorange  Ellipsoids  and  increasing  contribution  of  same  amount  of 


DOA  estimation  error  to  object  location  with  the  increase  in  range  sum . 4 

Eigure  4.  Baseline,  and  extended  baseline . 5 


Eigure  5.  Root  Sum  Square  (RSS)  of  range  sum,  baseline,  and  receive  angle  errors 
for  the  following  conditions;  Rt+Rr  =100km,  baseline  distance  (L)  =  10  km 


[3] . 6 

Eigure  6.  Bistatic  Geometry  and  Bistatic  Angle  (|3) . 15 

Eigure  7.  PCL  System . 19 

Eigure  8.  Minimum  Variance  Estimators  [4] . 26 

Eigure  9.  Phased  Array  Antenna  [17] . 31 

Eigure  10.  Direction  of  Arrival  Angles  with  respect  to  the  antenna . 35 

Eigure  11.  In  Phase  and  Quadrature  components  of  additive  noise  [16] . 36 

Eigure  12.  CRB  for  single  signal  case  when  SNR  =  OdB . 54 

Eigure  13.  The  CRB  as  of  a  moving  object  with  OdB  SNR  with  a  stationary  object 

of  the  same  SNR  located  at  0° . 57 

Eigure  14.  The  CRB  as  of  a  stationary  object  with  OdB  SNR  with  a  moving  object 
of  the  same  SNR . 59 


IX 


Figure  15.  The  CRBcm  of  the  moving  signal  source  with  OdB  SNR  with  a 
stationary  signal  source  of  the  same  SNR  when  the  Start  Phases  and  the 

Doppler  Shifts  of  the  two  signal  sources  are  the  same . 61 

Figure  16.  The  CRBcm  of  moving  signal  source  with  OdB  SNR  with  a  stationary 
signal  source  of  the  same  SNR  when  Start  Phases/ Doppler  Shifts  of  the 

moving,  and  stationary  signal  sources  are  [2,4]/ [8,-2]  respectively . 63 

Figure  17.  The  CRBcm  of  the  stationary  signal  source  with  OdB  SNR  with  a 
moving  signal  source  of  the  same  SNR  when  the  Start  Phases,  and  the 

Doppler  Shifts  of  the  two  signal  sources  are  the  same . 64 

Figure  18.  The  CRBcm  of  the  stationary  signal  source  with  OdB  SNR  with  a 

moving  signal  source  of  the  same  SNR  when  Start  Phases/ Doppler  Shifts  of 

the  moving,  and  stationary  signal  sources  are  [2,4]/ [8,-2]  respectively . 66 

Figure  19.  The  CRB  as  of  the  moving  signal  source  with  a  stationary  signal  source 
at  0°  and  a  transmitter  of  opportunity  at  70°  when  SNR  of  the  moving  and 

stationary  signal  sources  are  OdB  and  transmitter  SNR  50dB . 69 

Figure  20.  The  CRBcm  of  the  moving  signal  source  with  a  stationary  signal 

source  at  0°  and  a  transmitter  of  opportunity  at  70°  when  SNR  of  the  moving 
and  stationary  signal  sources  are  OdB,  and  transmitter  SNR  50dB;  and  Start 
Phases/ Doppler  Shifts  of  the  moving,  and  stationary  signal  sources, 

transmitter  are  [2,4,6]/ [8, -2,0]  respectively . 70 

Figure  21.  The  CRBcm  of  the  moving  signal  source  with  a  stationary  signal 

source  at  0°  and  a  transmitter  of  opportunity  at  12°  when  SNR  of  the  moving 


X 


and  stationary  signal  sources  are  OdB,  and  transmitter  SNR  50dB;  and  Start 
Phases/ Doppler  Shifts  of  the  moving,  and  stationary  signal  sources, 

transmitter  are  [2,4,6]/ [8,-2,0]  respectively . 71 

Figure  22.  The  CRB  as  of  the  Moving  Object  with  a  transmitter  of  opportunity  at 
different  locations  [-70°,  -50°,  -30°,  -15°,  0°]  when  SNR  of  the  moving  is  OdB, 

and  transmitter  SNR  50dB . 75 

Figure  23.  The  CRBcm  of  the  Moving  Object  with  a  transmitter  of  opportunity  at 
different  locations  [-70°,  -50°,  -30°,  -15°,  0°]  when  SNR  of  the  moving  is  OdB, 
and  transmitter  SNR  50dB;  and  Start  Phases/ Doppler  Shifts  of  the  moving, 
and  transmitter  are  [2,6]/ [8,0]  respectively . 78 


XI 


Notation 


Scalars,  Vectors,  Matrices 

Scalars  are  denoted  by  upper  or  lower  case  letters,  as  the  scalars  x  or  J. 

Vectors  are  denoted  by  underlined  letters,  as  the  vector  x  made  up  of 
components  Xi. 

Matrices  are  denoted  by  upper  case  letters  in  boldface  type,  as  the  matrix 
A  made  up  of  elements  Aij  (i*  row,  j*  column). 

E[.]  is  the  expectation  operator. 

.*  is  array  multiplication  operator. 

Superscripts 


Estimated,  Computed,  or  Measured  Value  (not  true  value). 
Matrix  Inverse. 

Vector  or  Matrix  Conjugate. 

Vector  or  Matrix  Transpose. 

Vector  or  Matrix  Hermitian. 

Real  part. 

Imaginary  Part. 


XU 


Acronym  List 


Abbreviation 

Acronym 

ACMA 

Algebraic  Constant  Modulus  Algorithm 

AFIT 

Air  Force  Institute  of  Technology 

BVR 

Beyond  Visual  Range 

CBF 

Conventional  Beam  Forming 

CM 

Constant  Modulus 

CPI 

Coherent  Processing  Interval 

CRB 

Cramer-Rao  Bound 

CW 

Continuous  Wave 

DOA 

Direction  of  Arrival 

DoD 

Department  of  Defense 

ESA 

Electronically  Steered  (Scanned)  Array 

FEBA 

Eorward  Edge  of  Battle  Area 

IEEE 

Institute  of  Electrical  and  Electronics  Engineers 

LOS 

Line  of  Sight 

ML 

Maximum  Likelihood 

MLE 

Maximum  Likelihood  Estimate 

MSE 

Mean  Squared  Error 

MUSIC 

Multiple  Signal  Characterization 

MV 

Minimum  Variance 

MVUB 

Minimum  Variance  Unbiased 

PCL 

Passive  Coherent  Location 

PDE 

Probability  Density  Eunction 

ULA 

Uniform  Linear  Array 

Xlll 


AFIT/GE/ENG/02M-23 


Abstract 

Direction  of  Arrival  (DOA)  estimation  of  signals  has  been  a  popular 
research  area  in  Signal  Processing.  DOA  estimation  also  has  a  significant  role  in 
the  object  location  process  of  Passive  Coherent  Location  (PCL)  systems.  PCL 
systems  have  been  in  open  literature  since  1986  and  their  applications  are  not  as 
clearly  understood  as  the  DOA  estimation  problem.  However,  they  are  the  focus 
of  many  current  research  efforts  and  show  much  promise. 

The  purpose  of  this  research  is  to  analyze  the  DOA  estimation  errors  in  a 
PCL  system.  The  performance  of  DOA  estimators  is  studied  using  the  Cramer- 
Rao  Bound  (CRB)  Theorem.  The  CRB  provides  a  lower  bound  on  the  variance  of 
unbiased  DOA  estimators.  Since  variance  is  a  desirable  property  for  measuring 
the  accuracy  of  an  estimator,  the  CRB  gives  a  good  indication  about  the 
performance  of  an  estimator. 

Previous  DOA  estimators  configured  with  array  antennas  used  the  array 
antenna  manifold,  or  the  properties  of  the  array  antenna  structure,  to  estimate 
signal  DOA.  Conventional  DOA  estimators  use  arbitrary  signal  (AS)  structures. 
Constant  Modulus  (CM)  DOA  estimators  restrict  the  input  signals  to  a  family  of 
constant  envelope  signals,  and  when  there  are  multiple  signals  in  the 
environment,  CM  DOA  estimators  are  able  to  separate  signals  from  each  other 
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using  the  CM  signal  property.  CM  estimators  then  estimate  the  DOA  for  each 
signal  individually. 

This  research  compares  the  CRB  for  AS  and  CM  DOA  estimators  for  a 
selected  system.  The  CRB  is  also  computed  for  this  system  when  single,  multiple 
and  moving  objects  are  present.  The  CRB  as  and  CRBcm  are  found  to  be  different 
for  the  multiple  signal  case  and  moving  object  cases. 
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STATISTICAL  ERROR  ANALYSIS  OF  A  DOA 
ESTIMATOR  FOR  A  PCL  SYSTEM  USING  THE 
CRAMER-RAO  BOUND  THEOREM 

1 .  INTRODUCTION 


1.1  Background 

The  leading  countries  in  history,  the  super  powers  of  their  own  time,  have 
always  showed  great  interest  in  scientific  research  and  development  and  its 
applications  in  military  matters.  The  lifetime  of  their  superiority  was  related  to 
how  long  they  were  able  to  maintain  their  position  as  scientific  leaders. 

Political  uncertainties  before  WWII  sent  a  message  of  the  upcoming 
conflict  and  every  nation  was  aware  that  there  were  going  to  be  many  firsts 
during  this  war,  like  the  first  use  of  the  Radio  Detection  and  Ranging  (RADAR) 
systems.  After  airpower  entered  the  warfare  arena,  early  detection  of  airborne 
objects  became  of  great  interest,  due  to  the  importance  of  quick  reaction  against 
enemy  airpower  [1].  In  the  last  century  there  have  been  numerous  research  and 
development  efforts  in  different  aspects  of  radar  systems,  causing  the  "lose  sight, 
lose  fight"  concept  to  go  Beyond  Visual  Range  (BVR)  of  human  eye.  Today,  there 
are  many  different  types  of  radar  systems  for  many  different  purposes  and  uses. 

Passive  Coherent  Location  (PCL)  and  the  related  terms  are  explained  in 
Chapter  2  but  basically  it  can  be  thought  as  a  passive  radar  system  which  is 
composed  of  only  a  receiver  unit.  Most  of  the  conventional  radar  concepts  still 
apply  to  PCL  because  the  illumination  of  the  object  still  takes  place;  not  by  a 
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transmitter  which  is  designed  for  this  purpose  but  by  a  transmitter  that  is 
already  in  that  environment.  Advantages  and  disadvantages  of  PCL  systems  are 
discussed  in  Chapter  2. 

1.2  Problem  Statement 

Monostatsic  radars  follow  the  generalizations  and  equations  that  are 
applicable  to  bistatic  radar.  Monostatic  object  location  errors  are  functions  of 
transmitted  pulse  width,  antenna  half  power  beam  width,  and  signal-to-noise 
ratio  [2]. 

There  are  inherent  errors  in  object  location  caused  by  bistatic  geometry  in 
addition  to  the  errors  that  are  caused  in  monostatic  radars.  The  most  basic 
system  architecture  scenario  of  a  PCL  system  is  the  bistatic  case,  which  is 
illustrated  in  Figure  1.  0t  is  the  transmit  angle  and  is  the  angle  from  the  reference 
direction  to  the  line  between  the  transmitter  and  the  object,  and  0r  is  the  receive 
angle  and  is  the  angle  from  the  reference  direction  to  the  line  between  the  receiver 
and  the  object.  Object  location  errors  that  apply  to  monostatic  radars  also  apply 
to  bistatic  radars  but  object  location  errors  in  bistatic  radars  have  additional 
dependence  on  the  receive  angle,  0r  [3]. 
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Figure  1.  Transmit  and  Receive  Angles  (0t,  0r). 

To  demonstrate  the  dependence  of  object  location  errors  on  the  receive 
angle  0r,  let's  introduce  a  basic  object  location  technique  for  a  bistatic  system 
from  [3].  Measuring  the  total  range  from  the  transmitter  to  object  (Rx)  and  from 
the  object  to  receiver  (Rr)  places  the  object  on  an  isorange  ellipsoid  as  shown  in 
Figure  2.  The  range  from  the  transmitter  to  an  arbitrary  point  on  this  isorange 
ellipsoid  plus  the  range  from  the  same  point  to  the  receiver  is  the  Range  Sum 
(Rt+Rr).  This  range  sum  is  constant  for  all  points  on  an  isorange  ellipsoid.  The 
foci  points  of  the  ellipsoid  are  the  locations  of  the  receiver  and  the  transmitter 
sites. 


T  arget 


Isorange  Ellipsoid 


Figure  2.  Isorange  Ellipsoid. 
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Given  range  sum  (Rt+Rr)  and  receive  angle  (0r),  Rr  is  calculated  by: 


Rr  = 


(R, +R,f-C 
2(R^  +  Rj^+ L  sin(0j^ )) 


(1) 


where  L  is  the  baseline  distance.  The  object  is  declared  to  be  at  range  Rr  and  its 


direction  is  determined  via  0r. 


Rx2  Rxi 


Figure  3.  Isorange  Ellipsoids  and  increasing  contribution  of  same  amount  of 
DOA  estimation  error  to  object  location  with  the  increase  in  range  sum. 


There  are  many  object  location  methods  in  PCL  systems,  with  each  based 
on  different  object  measurements.  The  most  basic  object  location  technique  is 
done  by  range  sum  measurement  and  DOA  estimation.  In  Figure  3  two  different 
isorange  ellipsoids  are  shown  for  two  different  range  sums.  DOA  estimation 
helps  us  determine  on  which  point  of  the  isorange  ellipsoid  the  object  is  located. 
Figure  3  illustrates  that  for  smaller  range  sums  the  error  in  DOA  estimation  has  a 
smaller  contribution  to  object  location  errors  and  as  the  object  range  increases  the 
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same  amount  of  angular  error  in  DOA  estimation  leads  to  an  increase  in  object 
location  errors.  When  the  estimated  object  direction  is  off  from  true  object 
direction  by  some  error  then  the  spatial  error  in  object  location  can  be 
represented  by  Rxi  and  Rx2,  assuming  that  range  sum  is  measured  correctly.  Rxi 
and  Rx2  in  Figure  3  are  cross-range  values  and  they  demonstrate  the  increasing 
contribution  of  same  amount  of  DOA  estimation  error  to  object  location  with  the 
increase  in  object  range  sum. 

In  Figure  5,  the  root  sum  square  errors  are  plotted  versus  receiver  look 
angle  (0r)  for  a  bistatic  system.  The  contribution  of  different  kinds  of  object 
location  errors  for  different  receive  angles  are  clearly  seen.  When  the  object  is  on 
the  extended  baseline  shown  in  Figure  4  (the  straight-line  that  passes  through  the 
transmitter  and  receiver,  excluding  the  line  segment  between  the  transmitter  and 
the  receiver),  transmit  and  receive  angles  are  equal.  This  means  that  when  an 
object  is  on  the  extended  baseline,  0t  which  is  equal  to  0r,  and  this  causes  the 
object  location  errors  to  converge  to  the  monostatic  case  [3]. 

Extended  ^ ^ ^  Baseline  ^  Extended 

Baseline  (L)  Baseline 

Transmitter  Receiver 

Figure  4.  Baseline,  and  extended  baseline. 
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RSS  of  Errors  vs.  Receive  Angle  Riot 


for  the  following  conditions:  Rt+Rr  =100km,  baseline  distance  (L)  =  10  km  [3]. 


The  receive  angle,  0r,  is  determined  by  Direction  of  Arrival  (DOA) 
estimation.  As  seen  in  Figure  5,  the  receive  angle  errors  have  the  greatest 
contribution  to  total  object  location  error  over  a  wide  range  of  receive  angles.  As 
the  receive  angle  approaches  +/-  90°,  the  contribution  of  receive  angle  error  to 
total  object  location  error  is  less  than  other  errors.  Section  1.5  introduces  the 
system  used  in  this  thesis  and  the  antenna  is  chosen  to  be  an  array  antenna.  Due 
to  array  antenna  characteristics,  reception  using  the  array  antenna  is  chosen  to  be 
in  the  approximate  region  from  -70°  to  70°.  This  region  is  where  receive  angle 
has  the  greatest  contribution  to  total  object  location  errors. 
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Keeping  in  mind  the  content  of  Figure  3  and  Figure  5,  clearly  any  increase 
in  DOA  estimation  accuracy  (receive  angle  estimation)  leads  to  more  accurate 
results  in  object  location.  The  curves  in  Figure  5  are  obtained  by  applying  the 
root  sum  square  (RSS)  concept  on  bistatic  object  location  technique.  It  is 
important  to  note  that  object  location  is  not  determined  in  this  fashion  for  all  PCL 
systems.  However,  these  figures  effectively  illustrate  the  importance  of  accuracy 
in  DOA  estimation  for  a  more  precise  estimate  of  object  location.  Since  DOA 
estimation  plays  a  major  role  in  object  location  of  PCL  systems,  the  accuracy  of 
DOA  estimation  is  the  primary  emphasis  of  this  thesis. 

1.3  Scope 

Error  analysis  of  PCL  systems  is  a  very  broad  research  area.  There  are 
electromagnetic  based  issues,  signal  processing  based  issues,  hardware  related 
issues,  software  related  issues,  etc.  This  thesis  focuses  on  the  signal  processing 
side  of  the  problem  and  performs  statistical  error  analysis.  As  mentioned  above, 
this  analysis  is  performed  on  the  most  essential  part  of  object  location,  the  DOA 
estimation. 

The  problem  is  considered  in  two-dimensions  (object  altitudes  are 
ignored).  When  the  problem  is  reduced  to  two-dimensions,  the  ellipsoids  in 
three-dimensional  cases  become  ellipses  in  two-dimensional  scenarios.  All 
objects  are  at  zero  altitude  and  the  DOA  estimation  problem  is  in  azimuth  only. 
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1.4  Approach/Methodology 

Different  statistical  approaches  in  spectral  estimation  have  introduced 
different  algorithms  in  signal  processing  for  estimating  signal  DOA.  The 
performance  of  these  algorithms  is  commonly  measured  using  the  variance  of 
the  DOA  estimate.  The  Cramer-Rao  Bound  (CRB)  is  a  statistical  theorem  that 
places  a  lower  bound  on  this  variance  for  any  unbiased  estimator,  given  the 
probability  density  function  of  the  error  is  the  same  among  those  estimators  [4, 
5], 

DOA  estimation  of  different  signals  using  an  array  antenna  has  become  a 
common  approach.  Most  algorithms  such  as  Multiple  Signal  Characterization 
(MUSIC),  Maximum  Likelihood  Estimation  (MLE)  [6]  and  Conventional  Beam 
Forming  (CBF)  use  the  properties  of  array  antennas  or  their  structure  and  have 
the  same  statistical  error  probability  distribution  function.  This  is  explained  in 
more  detail  in  Chapter  3. 

Improvements  in  spectral  estimation  enabled  the  use  of  signal  properties 
in  DOA  estimation,  in  addition  to,  the  array  manifold  or  properties  of  the  array 
structure.  In  [7]  an  algorithm  for  Constant  Modulus  (CM)  signals  has  been 
introduced  (signals  with  constant  envelope  using  phase  or  frequency 
modulation).  An  Algebraic  Constant  Modulus  Algorithm  (ACMA)  uses  the 
constant  modulus  property  of  signals  and  makes  a  source  separation  among 
signals  having  the  CM  property;  the  DOA  estimation  problem  is  decoupled  in 
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multiple  object  scenarios,  which  is  usually  the  case  in  real-world  applications. 
The  CM  DOA  estimators  have  a  different  CRB  and  this  is  explained  in  Chapter  3. 

This  thesis  discusses  the  derivation  and  calculation  of  the  CRB  for 
estimators  using  common  algorithms  over  a  wider  class  of  signals  (arbitrary 
signals)  and  the  CRB  for  estimators  using  CM  algorithms  (CM  signals). 
Therefore,  to  make  a  distinction  between  the  two,  the  CRB  for  Arbitrary  Signals 
(AS)  is  denoted  as  CRBas  and  the  CRB  for  Constant  Modulus  (CM)  signals  is 
denoted  as  CRBcm. 

The  difference  in  the  CRB  of  arbitrary  signals  and  CM  signals  is 
demonstrated  by  calculating  the  CRB  for  different  scenarios.  The  correlation  of 
incoming  signals  is  taken  into  account  throughout  the  calculations.  Changes  in 
the  CRB  with  respect  to  the  correlation  level  of  the  signals  is  shown  for  CRBas. 

1.5  System  Description 

The  CRB  is  directly  related  to  a  system  architecture;  thus,  a  basic  system 
model  is  established  for  calculating  the  CRB.  Since  all  the  algorithms  mentioned 
above  use  array  antennas,  the  model  consists  of  an  Uniform  Linear  Array  (ULA) 
antenna  with  16  elements  and  element  spacing  of  X/2.  In  each  Coherent 
Processing  Interval  (CPI),  512  realizations  (snap  shots)  are  considered. 
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1.6  Assumptions 

To  make  some  approximations  and  simplify  the  problem  at  hand,  this 
research  is  based  on  the  following  assumptions: 

•  This  research  assumes  the  interfering  noise  e(t)  is  circular  Gaussian 
distributed  white  noise  with  zero  mean  (E[e(t)]=0)  and  E[e(t)e^(t)]  =  al, 
where  e(t)  is  a  noise  vector  and  I  is  the  identity  matrix.  This  is  a 
reasonable  assumption  to  model  the  noise  in  the  environment. 

•  Objects  are  at  zero  altitude  and  DOA  estimation  is  performed  in  azimuth 
only.  This  assumption  simplifies  the  problem  being  faced  in  the  real 
world.  Therefore,  the  results  obtained  using  this  assumption  form  a 
baseline  to  extend  the  results  to  the  real  world  scenarios. 

•  Object  SNR  remains  constant  during  a  Coherent  Processing  Interval  (CPI). 
This  is  a  reasonable  assumption. 

•  Signal  DOA  remains  constant  during  a  CPI.  This  is  a  reasonable 
assumption  because  the  real-world  objects  can  be  maneuvering  objects. 
But  it  is  a  good  starting  point. 

•  The  antenna  is  calibrated  so  that  the  array  response  for  each  DOA  is 
known.  This  is  a  good  assumption  and  is  consistent  with  real-world 
applications. 

•  When  calculating  CRBcm,  all  signals  are  assumed  to  be  CM  signals. 
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•  Object  ranges  are  assumed  to  be  far  enough  for  the  signal  wave  fronts  to 
be  considered  as  plane  waves.  This  is  a  reasonable  assumption  since  this 
is  usually  the  case  in  radar  application  due  to  the  desired  object  detection 
ranges. 
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2. BACKGROUND 


The  purpose  of  this  chapter  is  to  explain  the  Passive  Coherent  Location 
(PCL)  system  and  its  related  issues.  The  basic  terms  related  to  PCL  are  explained 
thoroughly,  and  the  position  of  a  PCL  system  in  today's  Radar  world  is  made 
clear. 

Therefore,  this  chapter  first  discusses  developments  in  radar,  during 
which  basic  concepts  that  are  needed  for  a  good  understanding  of  a  PCL  system 
are  introduced.  Second,  the  PCL  system  is  explained,  including  PCL  system 
definition,  operation,  characteristics,  and  advantages  and  disadvantages  in  both 
general  and  military  applications.  Next,  the  Cramer  Rao  Lower  Bound  Theorem 
is  explained.  Finally,  phased  array  antennas  are  introduced  because  of  the  need 
to  establish  a  basic  understanding  of  the  "big  picture"  of  system  description  and 
system  operation. 

2.1  Development  of  RADAR  and  Basic  Terms 

When  Radio  Detection  and  Ranging  (RADAR)  systems  are  the  topic,  the 
detection  of  aircraft  or  airborne  objects  is  the  typical  consideration.  However, 
when  radar  was  first  developed  the  primary  objects  were  typically  naval  vessels. 
After  Hertz  showed  that  radio  waves  could  be  reflected  from  metallic  objects, 
and  formed  the  basis  for  the  concept  of  radar  [8],  Hulsmeyer  built  the  first 
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operational  radar  in  1904  to  detect  ships  [9].  Hulsmeyer's  radar  was  monostatic 
(collocated  receiving  and  transmitting  sites). 

Some  coincidences  also  contributed  to  the  development  of  radar.  In  1922, 
Taylor  and  Young  detected  a  change  in  the  signal  as  a  ship  passed  between  a 
transmitter  and  receiver  located  on  opposite  sides  of  a  river  [10].  This  example 
from  history  is  completely  compatible  with  today's  Continuous  Wave  (CW)  radar. 
CW  Radars  transmit  continuously  and  receive  while  transmitting  [Ij.  Although 
it  is  possible  to  build  a  CW  radar  with  a  collocated  receiver  and  transmitter,  CW 
radar  applications  are  much  easier  to  implement  with  a  bistatic  design  than  the 
monostatic  design. 

In  the  very  early  stages  of  radar  development,  detecting  ships  by  radar 
was  the  primary  concern  because  unlike  the  aircraft  of  that  time,  ships  were 
made  of  metal.  It  wasn't  until  the  late  1920's  that  aircraft  manufacturers  changed 
aircraft  designs  and  started  using  metal  on  the  fuselage  and  the  wings  (versus  a 
wooden  structure  with  fabric  redundant  cover)  [Ij.  After  this  change  in  aviation 
design,  detection  of  aircraft  by  radar  became  possible.  Each  country  began 
showing  interest  in  military  radar  applications  for  aircraft  detection  because  of 
the  political  imbalance,  which  lead  the  world  into  the  Second  World  War  (WW 
11).  The  revolution  in  technology  changed  the  tactics  of  war  that  would  be 
fought.  Meanwhile,  aviation  as  a  military  power  became  of  primary  importance. 
These  reasons  caused  radar  advancement  to  exceed  a  natural  development  line 
and  drew  a  great  deal  of  attention  to  this  field. 
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In  bistatic  radars,  the  transmitter  and  the  receiver  are  separated,  but  the 
distance  between  the  transmitter  and  the  receiver  is  known  (Figure  6).  The 
distance  between  the  transmitter  and  the  receiver  sites  is  called  the  baseline 
distance  (L);  requirements  on  the  length  of  the  baseline  distance  are  vague  in  most 
texts.  The  IEEE  description  of  bistatic  radar  also  does  not  completely  describe 
this  baseline  distance  [3].  Skolnik  provides  better  explanation;  he  argues  that  the 
distance  should  change  relative  to  the  object  distance.  A  baseline  can  be 
relatively  small  for  objects  that  are  close  in  range.  On  the  other  hand  it  must  be 
larger  for  the  long-range  objects.  A  different  way  of  describing  this  baseline 
requirement  is  the  need  for  two  different  electromagnetic  paths  for  the 
propagation  of  the  transmitted  signal  from  transmitter  to  the  object  and 
propagation  of  the  object  echo  from  the  object  to  the  receiver  [11]. 

Bistatic  angle  (|3)  is  equal  to  the  difference  between  the  transmit  and  the 
receive  angles  0t  and  0r.  0t  and  0Rare  defined  in  Section  1.2. 

P  =  0T-0R  (2) 

Monostatic  radar  can  be  considered  as  a  special  case  of  bistatic  radar  with 
Rt  equal  to  Rr,  0t  equal  to  0r,  and  L  =  0.  Since  0t  equals  0r,  |3  is  equal  to  zero  in 
monostatic  radar.  In  bistatic  radar  applications,  the  bistatic  angle  becomes  zero 
as  the  object  approaches  the  extended  baseline;  this  is  called  over  the  shoulder 
operation. 
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Figure  6.  Bistatic  Geometry  and  Bistatic  Angle  (p). 


Most  of  the  earliest  radars  were  bistatic.  In  some  early  applications, 
systems  were  built  to  detect  objects  on  the  baseline  between  the  transmitter  and 
the  receiver.  This  special  case  is  called  the  forward  scatter  case  (or  forward  scatter 
geometry).  In  a  forward  scatter  case,  the  object  echo,  which  is  scattered  forward 
from  the  object  towards  the  receiver,  is  of  interest  rather  than  the  object  echo  that 
is  reflected  back  towards  the  transmitter.  Initially,  bistatic  radars  that  could 
detect  objects  crossing  the  baseline  using  the  forward  scatter  from  the  objects 
were  developed,  and  they  were  called  forward  scatter  fences.  Japan,  the  Soviet 
Union,  and  France  had  approximately  200  of  these  forward  scatter  fences  before 
and  during  WW  II  [3].  After  forward  scatter  fences,  use  of  bistatic  radar  was 
extended  out  of  the  baseline  area.  The  Klein  Heidelberg  device,  the  first  bistatic- 
of-opportunity  system,  was  used  by  Germans  in  WW  II  [12]. 
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Object  detection  and  location  is  a  much  more  complicated  task  in  bistatic 
radars,  as  compared  to  monostatic  radars.  Thus,  upon  the  invention  of  the 
duplexer,  a  device  enabling  the  use  of  a  single  antenna  for  both  transmission  and 
reception  by  avoiding  transmitter  and  receiver  interference,  monostatic  radars 
gained  greater  popularity  over  bistatic  radars,  especially  after  WW  II  [13].  This 
lead  to  a  drastic  reduction  in  bistatic  radar  research. 

Bistatic  radar  establishes  a  basic  understanding  for  multistatic  radars. 
When  more  than  one  receiver  or  more  than  one  transmitter  is  used  for  object 
detection,  and  object  detection  is  done  coherently  by  establishing  a  network 
between  the  receiver  sites,  the  system  established  is  called  a  multistatic  system  [2]. 
Each  transmitter  and  receiver  pair  of  a  multistatic  system  is  a  bistatic  system. 
Creating  a  network  of  these  basic  bistatic  elements  turns  these  separate  bistatic 
systems  into  a  multistatic  system. 

Each  radar  measurement  from  a  different  transmitter-receiver  pair  (or 
basic  bistatic  element)  is  sent  to  a  central  station  where  data  processing  is 
performed.  There  is  a  natural  trade-off  between  the  coverage  area  of  the 
multistatic  system  and  the  object  estimates.  Creating  a  larger  common  coverage 
area  of  the  bistatic  elements  of  a  multistatic  system  can  enhance  object  estimates, 
but  this  reduces  the  total  coverage  area  of  the  multistatic  system.  If  a  wider 
coverage  area  is  desired,  then  the  basic  bistatic  elements  must  be  spread  apart, 
thus  reducing  the  common  coverage  area.  This  in  turn  degrades  the  object 
estimation  [3]. 
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2.2  Passive  Coherent  Location  System 

There  was  a  turn-around  on  the  reduction  of  bistatic  radar  research 
around  1980s.  Some  applications  of  bistatic  radar  became  really  popular,  e.g., 
semi-active  missiles.  Improvements  in  destroying  enemy  radars  lead  to  the 
development  of  anti-radiation  missiles  (ARM).  Bistatic  radars  are  less  susceptible 
to  ARMs  by  locating  the  receiver  close  to  the  Forward  Edge  of  Battle  Area 
(FEB A)  and  the  transmitter  further  away  in  a  more  secure  environment,  well  out 
of  ARMs  reach.  Basically  in  some  scenarios  bistatic  radars  have  advantages  over 
monostatic  radars.  This  realization  made  the  bistatic  radars  an  interesting 
research  topic  once  again. 

Bistatic  radars  use  regular  radar  transmitters  designed  for  radar 
applications.  As  the  researches  on  bistatic  radars  kept  going  new  ideas  such  as 
using  different  transmitters  were  considered  and  that  is  how  idea  of  Passive 
Coherent  Location  (PCL)  was  born. 

PCL  systems  use  ambient  radio  waves  as  a  source  of  illumination  in  order 
to  detect  and  track  objects.  The  most  common  type  of  radio  waves,  available  all 
around  the  world  for  this  purpose,  are  television  and  PM  radio  broadcasts,  even 
though  they  might  differ  somewhat  from  one  country  to  another.  A  PCL  system 
can  be  either  multistatic  or  bistatic,  depending  on  the  kind  of  network  and 
number  of  transmitters  and  receivers  used  in  the  system.  The  most  basic  form  of 
a  PCL  system  is  the  bistatic  case  when  only  one  transmitter  and  receiver  are  used 
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for  object  detection  and  tracking  [2],  Figure  7  shows  the  basic  architecture  of  a 
PCL  system. 

2.2.1  Types  o£  PCL 

There  are  three  types  of  PCL  systems,  classified  with  respect  to  the  kinds 
of  signal  that  they  utilize  to  obtain  object  estimations.  These  are  Narrow  band 
PCL,  Wide  band  PCL,  and  Pulsed  PCL  [14]. 

Narrow  band  PCL  uses  the  narrow  video  carrier  or  audio  carrier  portion  of 
a  TV  waveform  where  most  of  the  signal  energy  is  contained.  This  type  of  PCL 
enables  Doppler  and/ or  DO  A  measurements. 

Wide  band  PCL  uses  the  modulation  spectrum  of  an  FM  waveform,  which 
has  a  broader  band.  The  methods  used  in  Wide  band  PCL  enable  range  and/ or 
DOA  measurements.  In  some  cases  Doppler  measurements  are  also  obtained. 

Pulsed  PCL  system  uses  pulsed  waveforms  from  pulsed  radars. 

2.2.2  System  Operation 

Most  PCL  systems  receive  a  direct  signal  from  the  transmitter  and  another 
signal  scattered  from  the  object.  Thus  the  correlation  between  the  scattered 
signal  from  the  object  and  direct  (reference)  signal  will  give  the  object's  range 
estimate.  If  the  reference  signal  is  distorted  by  any  means,  this  will  present  some 
deviation  from  the  original  object  information.  As  Howland  states  in  [2], 
"Accurate  bearing  information  can  be  obtained  using  two  or  more  antennas, 
using  techniques  such  as  phase  interferometry,  amplitude  monopulse, 
conventional  beam  forming,  and  super  resolution." 
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There  are  times  when  the  direct  signal  from  the  transmitter  interferes  with 
the  beam  looking  in  the  object  direction  and  the  receiving  object  information. 
When  this  occurs  one  of  the  simplest  methods  to  get  rid  of  this  direct  signal 
interference  is  to  use  the  reference  signal  to  cancel  it  out.  The  Adaptive  beam 
forming  technique  can  be  used  to  null  out  most  interference  that  might  be 
encountered  [15]. 


Figure  7.  PCL  System. 

2.2.3  Characteristics  of  PCL 

PCL  is  called  a  "coherent"  system  because  it  uses  both  the  phase  and  the 
amplitude  of  the  received  waveform  to  obtain  object  estimates.  PCL  measures 
range,  Doppler,  and  DOA  parameters  of  the  object  when  bistatic  architecture  is 
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used.  In  Multistatic  cases.  Time  Difference  of  Arrival  (TDOA)  and  Differential 
Doppler  measurements  are  also  of  interest  [2], 

A  PCL  system  uses  a  non-cooperative  transmitter,  which  means  the 
system  has  no  control  of  the  transmitter.  Thus,  these  transmitters  are  also  called 
transmitters  of  opportunity.  A  system  can  be  designed  to  use  different 
waveforms,  but  it  is  restricted  to  use  whichever  waveform  exists  in  its 
environment. 

Because  PCL  is  bistatic,  the  transmitter  and  the  receiver  must  both  have 
clear  line  of  sight  (LOS)  to  the  object.  This  requirement  has  to  be  taken  into 
consideration  for  the  system  deployment.  Because  there  is  no  control  on  the 
transmitter  of  opportunity,  a  receiver  must  be  placed  to  satisfy  this  requirement. 

Since  PCL  uses  TV  and  FM  radio  broadcasts,  which  are  not  designed  for 
radar  processing  purposes,  these  waveforms  have  poor  radar  ambiguity  function 
[14].  For  that  reason,  a  PCL  system  must  use  a  non-linear  estimation  or 
optimization  algorithm  to  resolve  object  range,  bearing,  heading  and  velocity,  as 
opposed  to  conventional  radars  [16]. 

A  PCL  system  differs  from  a  conventional  monostatic  radar  system  by  the 
frequency  spectrum  it  uses,  and  its  bistatic  or  multistatic  geometry  [14]. 
Multistatic  geometry  enables  additional  object  data,  and  as  a  result  of  this  much 
more  accurate  object  location  estimates  can  be  performed.  A  PCL  system  must 
operate  in  VHF  and/  or  UHF  frequency  range  in  units  of  MHz;  this  is  the 
frequency  range  used  for  all  TV  and  FM  radio  transmissions.  Most  conventional 
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radars  operate  in  the  GHz  range.  Therefore,  processing  techniques  that  are  used 
with  today's  conventional  radar  must  be  modified;  this  is  an  arduous  task,  which 
can  be  done  only  by  the  extensive  use  of  digital  processing. 

2.2.4  Advantages  of  PCL 

PCL  is  potentially  a  low  cost  system.  It  uses  transmitters  of  opportunity 
as  sources  of  illumination;  the  radiation  sources  are  already  out  there  in  the 
environment.  This  saves  PCL  users  the  money  that  will  be  spent  on  the 
transmitter  [16].  On  the  other  hand,  the  receiver  must  be  more  complex  than 
conventional  radar  to  compensate  for  operating  on  a  "free"  transmitter,  but  all  of 
the  modifications  can  be  done  in  software.  The  hardware  requirements  for  the 
PCL  system  are  the  same  as  the  equipment  contained  in  the  receiver  unit  of 
conventional  radar;  however,  this  equipment  has  to  be  designed  to  operate  in  the 
frequency  spectrum  that  PCL  uses. 

Therefore,  emphasis  is  on  the  software  in  PCL  systems.  Software  does  the 
correlation  between  the  reference  signal  and  the  scattered  signal  to  estimate  the 
range,  the  Doppler  processing,  and  the  digital  beam  forming.  The  future 
improvements  of  the  system  seem  to  be  in  more  robust  software  rather  than  new 
hardware  requirements.  The  possible  mid-life  upgrade  on  a  PCL  system  will 
most  probably  be  cost  effective. 

The  frequency  separation  among  PCL  systems  is  not  as  big  an  issue  as  it  is 
in  conventional  radars  [2].  Conventional  radars  are  designed  to  transmit  and 
receive  at  one  or  more  frequencies.  Today  there  are  numerous  radar  designs. 
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and  most  of  these  radars  are  likely  to  have  the  same  operating  frequencies,  thus 
causing  an  interference  problem  among  radars  using  similar  frequency  ranges. 
This  is  not  the  case  in  PCL.  VHP  and  UHF  frequencies  have  many  users  with 
frequency  allocation  for  each  user;  more  than  one  PCL  receiver  can  operate  out  of 
the  same  transmitter  in  multistatic  cases.  The  interference  problems  in  PCL 
occur  by  the  reception  of  a  different  transmitter  or  scattering  of  different 
transmitter  broadcasts  from  different  objects  operating  at  the  same  frequency. 
Once  a  PCL  system  is  deployed  to  an  area,  and  its  set-up  is  complete,  these  kinds 
of  interferences  can  be  experienced.  Appropriate  receiver  placement  to  exploit 
the  terrain  can  be  used  to  prevent  interfering  transmitter  broadcasts,  but  object 
returns  from  these  transmitter  signals  are  harder  to  eliminate.  These  returns 
must  be  compensated  for  in  the  software. 

2.2.5  Disadvantages  o£  PCL 

The  most  essential  disadvantage  of  a  PCL  system  is  having  no  control 
over  the  transmitting  site.  This  can  be  considered  in  two  ways.  The  first  is 
having  no  physical  control  over  the  transmitter:  when  to  transmit,  etc.  The 
second  is  having  no  control  on  the  transmission:  the  broadcast  frequency, 
average  power,  the  transmitted  waveform,  etc.  Since  commercial  broadcasts  are 
not  designed  for  radar  processing,  they  have  a  poor  ambiguity  function. 

Other  disadvantages  are  caused  by  the  bistatic  nature  of  the  system.  In 
order  to  detect  an  object  in  a  surveillance  region  both  the  transmitter  and  the 
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receiver  must  have  line  of  site  to  the  desired  surveillance  region.  When  terrain  is 
a  factor  this  might  cause  severe  limitations  for  low  altitude  objects. 

There  are  also  some  limitations  on  the  PCL  antenna.  Small  antenna  size  is 
preferred  due  to  the  efforts  of  making  the  system  low  cost,  mobile,  and  difficult 
to  detect.  The  trade  off  for  a  smaller  antenna  size  is  reduced  detection  range, 
besides  a  degraded  spatial  resolution. 

Some  PCL  receivers  have  very  high  receiver  noise  figures.  This  reduces 
thermal  noise-limited  detection  ranges  usually  by  a  factor  of  two,  or  even  more. 
Interference  of  Galactic  noise  is  one  of  the  reasons  for  this  reduction.  In  [1], 
Skolnik  uses  the  term  cosmic  noise  instead  of  Galactic  noise,  and  describes  it  as  a 
continuous  background  caused  by  electromagnetic  radiations  from 
extraterrestrial  sources  in  our  own  galaxy,  extragalactic  sources,  and  radio  stars. 
Galactic  noise  decreases  with  increasing  frequency.  Hence,  it  is  negligible  at 
frequencies  above  UHF.  However,  the  frequency  range  of  most  PCL  systems  is 
still  affected  by  the  Galactic  noise;  this  contributes  to  the  increased  noise  figure  of 
the  PCL  receivers.  Having  other  transmitters,  or  the  original  transmitter  of 
opportunity  in  the  surveillance  region  can  also  contribute  to  the  higher  noise 
figure  in  PCL  receivers  [14]. 

After  giving  the  definition  to  PCL,  mentioning  its  characteristics, 
advantages,  and  disadvantages  now  it  is  time  to  touch  on  the  CRB  theorem 
which  is  used  in  this  thesis  for  characterizing  the  general  performance  of  the 
DOA  estimators  used  for  PCL  purposes. 
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2.3  Cramer-Rao  Lower  Bound  (CRLB,  CRB) 

The  analysis  of  DOA  estimation  uses  the  statistical  technique  Cramer-Rao 
Bound  (CRB),  which  determines  the  best  achievable  estimate  using  the  lower 
bound  on  the  error  variance.  Therefore,  the  statistics  of  the  measurement  errors 
must  be  known  or  modeled. 

In  this  research  the  Fisher  information  matrix  depends  on  the  likelihood 
function,  which  is  expressed  in  terms  of  variance  of  errors  in  the  each  phase 
measurement  of  all  the  array  antenna  elements,  o.  Therefore,  when  the  Fisher 
matrix  is  inverted,  the  result  of  the  CRB  is  a  function  of  the  measurement  error 
variance,  o. 

For  this  reason,  CRB  enables  the  evaluation  criteria  for  the  marginal 
measurement  errors.  From  that  point,  accuracy  of  an  estimate  of  an  object's 
direction  can  be  determined  depending  upon  measurement  errors  and  the 
original  direction  of  the  object. 

When  a  parameter  is  being  estimated,  the  properties  of  the  estimator  - 
such  as  variance,  or  bias-  are  usually  of  great  concern.  Greater  insight  for  the 
estimator  provides  better  insight  about  the  actual  estimate  of  the  parameter  of 
interest,  in  particular  the  reliability  and  accuracy  of  the  estimate.  More  detailed 
discussion  on  these  terms  is  available  in  [4,  5]. 
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2.3.1  Unbiased  Estimators 

If  the  expected  value  of  the  estimator  output  is  the  true  value  of  the 
parameter  of  interest  then  this  estimator  is  unbiased.  Mathematically  an  estimator 

0  is  unbiased  if; 

£[01=0,  a<0<b  (3) 

where  0  is  the  unknown  parameter,  (a,b)  is  the  range  of  possible  values  for  0,  0  is 

an  estimator  for  the  parameter  0,  and  E[.]  is  the  expected  value  operator. 

2.3.2  Minimum  Variance  Unbiased  (MVUB)  Estimator 

Trying  to  design  or  select  an  estimator  from  a  family  of  estimators  is 
difficult.  Because  of  this  some  criteria  are  needed  to  compare  how  optimal  each 

estimator  is.  Mean  Square  Error  (MSE)  between  the  estimate  0  and  the  true  value 
0  is  one  of  them  and  it  is  defined  as: 

mse  (§)  =  E  (  0  -  0  )  (4) 

MSE  can  also  be  rewritten  as: 

mse  (§)  =  E  I  (  0-e[  §]  )  +(  e[  §]  -0  )  "|  (5) 

mse  (§)  =  e|  I  0  -E  0  j  |  +2{e  0  -  gg)(E  0  -  oj  +  -  O) 

(6) 

mse  (§)  =  var^oj  +  bias^^o)  (7) 
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where  Ha  =  E 


Since  |a,-  -E  0  =0 ,  the  whole  second  term  in  Equation  (6)  goes  to  zero  and 


the  remaining  part  can  be  written  as  a  function  of  the  variance  and  the  bias  of  the 
estimator. 


In  an  unbiased  estimator,  the  bias2(0)  term  in  Equation  (7)  is  zero.  In  this 

case,  MSE  is  only  a  function  of  the  variance  and  is  minimized  by  minimizing  the 
variance  for  the  unbiased  estimators. 

In  order  to  test  whether  an  estimator  is  unbiased  or  not,  its  expected  value 
is  checked.  If  the  expected  value  of  an  estimator  is  equal  to  the  unknown 
parameter  0,  then  this  estimator  is  unbiased  as  in  Section  2.3.1.  The  discussion  of 
Minimum  Variance  (MV)  is  restricted  to  the  class  of  unbiased  estimators. 


Figures.  Minimum  Variance  Estimators  [4]. 

Figure  8a  shows  the  variance  of  three  different  unbiased  estimators.  Since 

estimator  03  has  the  lowest  variance  for  all  0  in  Figure  8a,  it  is  the  MV  estimator. 


On  the  other  hand  in  Figure  8b  03  is  the  MV  estimator  when  0<a  and  0>b,  and  02 
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is  the  MV  estimator  when  a  <  0  <  b.  In  this  case,  there  is  no  single  MV  estimator 
for  all  0. 

A  Minimum  Variance  Unbiased  (MVUB)  estimator  is  an  unbiased  estimator 
that  has  the  minimum  variance  for  all  the  unknown  parameter  0.  Sometimes  no 
MVUB  estimator  can  be  found  in  the  search  for  an  estimator.  MVUB  estimators 
do  not  always  exist.  The  MSB  of  an  MVUB  estimator  is  equal  to  its  variance. 

2.3.3  Likelihood,  Log-likelihood,  and  Score  Functions 

Observed  data  enables  the  estimation  of  the  unknown  parameter.  Since 
observed  data  is  dependent  on  the  pdf  of  the  data,  estimation  starts  with  the 
determination  of  the  pdf  for  that  data.  If  the  pdf  of  the  observed  data  is 
independent  of  the  desired  parameter,  then  that  parameter  cannot  be  estimated 
using  that  set  of  data.  The  pdf  of  a  random  variable  X  parameterized  by  the 
unknown  parameter  0  is  denoted  as  f0(X). 

A  simple  example  makes  it  easier  to  understand  this  concept.  Let  x  be  the 
observed  or  measured  data,  n  be  additive  white  gaussian  noise  with  zero  mean 
and  variance  o^,  and  0  be  the  unknown  parameter.  The  relation  between  these 
terms  are  formulized  as 

x[t]  =  y  +  n[t]. 

In  this  case  a  good  estimator  for  0  is  0  =  x[t].  Since  E{0}  =  E{x[t]}  =  0,  0  is  an 


unbiased  estimator.  Its  variance  is  and  the  accuracy  of  the  estimator  gets 


better  as  decreases  and  the  pdf  is  f  (x)= — =J= 


exp 


a 


2o 
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"When  the  pdf  is  viewed  as  a  function  of  the  unknown  parameter  (with  x  - 
the  data-  fixed)  it  is  termed  the  Likelihood  Function"  [4],  Since  X  is  a  random 
variable  the  change  in  X  is  random.  Fixing  X  is  done  via  inputting  a  set  of 

observation  data  to  the  joint  pdf  of  X  and  0  (f0(X)).  x  is  denoted  as  the  observed 

data  and  the  mathematical  representation  of  the  likelihood  function  is: 

i(e,x)  =  f.  (ic)  (8) 

and  the  Log-Likelihood  Function  is  simply  the  natural  log  of  the  likelihood 
function: 

l(0,x)  =ln[l(0,x)]  =ln[  (x)  ].  (9) 

Differentiating  the  log-likelihood  function  with  respect  to  0  gives  the  score 
function  [5],  whose  mathematical  representation  is 

‘‘(e-x)  =  f  L(e,x)  =  ^in[i(e,x)]  (lO) 

The  maximum  likelihood  estimate  (MLE)  of  0  is  found  by  setting  the  score 
function  equal  to  zero  and  solving  for  0. 

s(0,x)=0  (11) 

This  is  equivalent  to  finding  the  maximum  of  fg  (xj . 

Using  X  instead  of  x,  which  represents  the  observed  data  is  usually  done 
for  notational  convenience  and  this  notation  is  also  used  in  this  thesis. 
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2.3.4  Cramer-Rao  Bound  Theorem 


The  Cramer-Rao  Bound  is  a  lower  bound  on  the  covariance  matrix  of  any 
unbiased  estimator  of  parameter  0.  Thus,  no  unbiased  estimator  can  exist  whose 
variance  is  lower  than  the  CRB  for  each  value  that  0  can  take  in  the  parameter 
space.  An  unbiased  estimator  that  attains  the  CRB  is  by  definition  the  MVUB 
estimator.  Thus,  the  CRB  is  a  great  aid  to  determine  if  an  unbiased  estimator  is 
the  MVUB.  It  sets  a  benchmark  against  which  we  can  compare  the  performance 
of  an  unbiased  estimator. 

As  mentioned  above,  values  of  0  that  makes  the  score  function  equal  to 


zero  are  Maximum  Likelihood  (ML)  estimates  for  0.  The  score  function  has  mean 
zero;  in  [4]  this  is  called  the  Regularity  Condition  and  mathematically: 


E[s(0,X)]  =  E 


T^lnfe(X) 

30  ^ 


=  J  1  f.  4‘ = “  (12) 
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The  Fisher  Information  Matrix  (  J(0)  )  is  the  covariance  matrix  of  the  score 
function: 


J(0)  =  e[s(0,X)s^(0,X)]  =  E 


(13) 


Some  simplifications  in  (13)  lead  to  another  formula  for  the  Fisher  Information 
matrix  [5]: 
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J(0)  =  -E 


d  f  d  ^ 

—  —  lnf„(X) 

aeUe  ) 


(14) 


The  individual  elements  of  the  Fisher  Information  matrix  can  be  found  by: 


Ji.=E 


Ainf  (x)^lnf9(X) 


Infe(X) 


ae^aBj 


(15) 


The  Cramer-Rao  Bound  theorem  indicates  that  each  diagonal  element  of 

the  covariance  matrix  (Ce)  of  an  unbiased  estimator  0  is  greater  than  or  equal  to 

the  corresponding  element  of  the  matrix  obtained  by  inversing  the  Fisher 
Information  matrix.  This  also  means  that  J  (0)  has  to  be  invertible.  That  is, 

q  2  J-'  (e)  (16) 


The  (i,i)  element  of  Cg  is  the  mean-squared  error  of  the  0i  estimate  of  the 
unknown  parameter  0i  and  the  (i,i)  element  of  J"^(0)  lower  bounds  it: 


C.  =E 


(e-et]>j:‘(e). 


(17) 


An  efficient  estimator  is  an  unbiased  estimator  whose  covariance  matrix  is 
equal  to  the  inverse  of  Fisher  Information  matrix  [5]: 

Cj  =  J"  (e) .  (18) 

Fach  efficient  estimator  is  the  MVUB  estimator  but  the  opposite  is  not 

true. 
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2.4  Phased  Array  Antennas 

Phased  Array  Antennas  are  mounted  on  a  steady  platform  and  their  beam 
is  steered  electronically  instead  of  moving  the  platform  of  the  antenna  to 
mechanically  change  the  direction  of  the  beam.  Controlling  the  phase  difference 
of  the  radiating  elements  both  at  transmission  and  reception  does  the  electronic 
steering  (Figure  9).  Therefore  in  some  literature  these  antennas  are  also  called 
Electronically  Steered  Array  Antennas  (ESAs). 

[1,  17,  18,  19]  present  more  detailed  information  on  phased  array 
antennas;  here  the  relation  between  the  phase  delay  between  the  radiating 
elements  and  the  desired  angle  of  transmission/  reception  is  shown. 


Figure  9.  Phased  Array  Antenna  [17]. 

In  Figure  9,  it  can  be  seen  that  the  radiating  element  in  point  A  must  lead 
the  other  radiating  elements.  Each  element  leads  the  one  that  is  located  on  its 
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right  side  (because  of  the  way  it  is  illustrated  in  Figure  9).  When  all  the  elements 
of  an  array  antenna  are  uniformly  separated,  this  type  of  antenna  is  called 
Uniform  Linear  Array  (ULA).  In  a  ULA  the  phase  difference  between  any  of  the 
two  elements  located  next  to  each  other  is  the  same. 

The  total  change  in  phase  over  a  distance  of  one  wavelength  (X)  is  271.  The 
ratio  of  2ti  to  X  gives  the  phase  per  unit  distance.  Phase  difference  ( AO )  over  a 
specific  distance  (AR)  is: 

AO  =  —  AR  (19) 

As  can  be  seen  from  Figure  9,  AR  is  equal  to  AR  =  (n-l)xdsin(a).  Substituting 
this  into  (19)  yields: 

271 

AO  =  — X  (n  -  l)d  sin(a)  (20) 

where  a  is  the  angle  off  broadside  of  the  antenna,  the  angle  between  the  antenna 
broadside  direction  and  the  direction  that  the  main  beam  is  pointing  at. 

AO  is  the  total  phase  difference  between  the  first  and  last  elements. 
Element  to  element  phase  difference  is  the  same  between  all  the  elements  of  a 
ULA  and  it  can  be  calculated  by  (20)  by  substituting  2  for  n.  Result  of  this  is  the 
phase  difference  between  the  first  and  the  second  elements,  which  is  equal  to  the 
phase  differences  between  all  the  adjacent  elements. 

The  phase  differences  calculated  using  the  formulas  above  hold  in  theory 
but  in  real  life  this  might  not  be  the  case.  The  near  field  interaction  of  the 
radiating  elements  (also  called  mutual  coupling  of  the  elements)  affects  the 
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radiation  pattern.  Therefore  each  antenna  must  be  calibrated  in  order  to  steer  the 
beam  to  the  exact  angle  off  the  antenna  broadside  direction. 

2.5  Summary 

This  chapter  explained  the  terms  required  for  a  better  understanding  of 
PCL  systems  and  discussed  PCL  characteristics,  advantages,  and  disadvantages. 
The  CRB  theorem  was  discussed,  and  the  use  of  array  antennas  in  the  DOA 
estimation  applications  was  shown.  All  of  these  sections  mentioned  above 
provide  insight  for  the  methodology  outlined  in  Chapter  3. 
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3 .  METHODOLOGY 


As  mentioned  in  Section  2.4  the  phase  differences  between  the  elements 
under  noise  free  conditions  results  with  Equation  (20).  But  the  noise  measured 
along  with  the  signals  affects  the  measured  phase  differences  between  the  array 
elements  and  therefore,  additive  noise  contributing  to  the  measurements  causes 
the  error  in  DO  A  estimation  when  using  array  antennas.  In  order  to  determine 
the  CRB  for  an  estimator,  the  statistical  distribution  of  the  error  or  pdf  of  the 
interfering  noise  must  be  known.  If  these  statistics  are  unknown,  then  they  must 
be  modeled  as  accurately  as  possible. 

3.1  System  Model 

As  mentioned  shortly  in  Chapter  1  a  Uniform  Linear  Array  (ULA)  antenna 
is  modeled  whose  element  spacing  is  half  a  wavelength  or  X/2.  The  number  of 
antenna  elements  is  m  -16;  the  number  of  realizations  (or  snap  shots)  N  =  512.  n 
is  the  number  of  signals  impinging  on  the  antenna. 

Figure  10  illustrates  how  the  DO  A  angles  are  referenced  throughout  this 
thesis.  0°  represent  the  antenna  broadside  direction;  the  angle  is  measured  in  a 
clockwise  fashion,  90°represents  a  one  quarter  rotation  to  the  right,  while  -90°  is 
a  rotation  to  the  left.  Thus,  -90°  and  90°  can  also  be  referred  as  90°  left  and  90° 
right. 
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DOA  of  0° 


▲ 


DOA  of  ^ -  -  --  --  --  --  --  --  --  -  - DOA  of 

-90°  V _ _ ^  90° 


16  Element  ULA 

Figure  10.  Direction  of  Arrival  Angles  with  respect  to  the  antenna. 

3.2  Derivation  of  CRB  for  Arbitrary  Signals  (CRB as) 

Derivation  of  CRB  as  for  DOA  follows  [6] .  Measured  data  can  be  modeled 
by, 

y(t)=  A(0R)x(t)+e(t)  t  =  l,2,---,N  (21) 

where 

y(t)  is  a  m-by-1  data  vector  including  measurement  noise, 
x(t)  is  a  n-by-1  vector  and  its  elements  are  signal  amplitudes, 
e(t)  is  the  m-by-1  noise  vector, 

A  is  a  m-by-n  array  response  matrix,  which  is  given  by  Equation  (22). 

A(0j  =  [a(0,j  -  a(0,J]  (22) 

Each  a(0Ri)  in  Equation  (22)  is  a  steering  vector  of  the  i*  signal  and  contains  the 
phase  delays  of  the  corresponding  elements;  0r  is  a  n-element  receive  angle 
(DOA)  vector  whose  elements  are  real  parameters  [0ri,  ...  ,  0Rn]^.  One  must  be 
cautious  with  notation.  0  all  by  itself  represents  the  unknown  parameters,  as 
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mentioned  in  the  definition  of  CRB  in  Section  2.3,  and  0r  represents  the  DOA  or 
receive  angle. 

Figure  11  shows  how  additive  noise  affects  the  measured  signal  in  one 
channel  of  the  array.  The  in-coming  signal  is  corrupted  by  an  in-phase  (real  part, 

ei(t))  and  quadrature  (imaginary  part,  ei(t))  component  of  the  noise  (ei(t)  indicates 


the  i*  element  of  noise  vector  e(t)).  The  resultant  phase  error  is  Vj/e  which  causes 
DOA  estimates  to  deviate  from  the  true  value. 


Quadrature  phase 
component  of 
noise 


Figure  11.  In  Phase  and  Quadrature  components  of  additive  noise  [16]. 


From  Equation  (21)  e(t)  can  be  written  as  in  Equation  (23): 

e(t)=y(t)-Ax(t)  (23) 

In  Section  1.6,  each  noise  vector  e(t)  is  assumed  to  have  a  circular  Gaussian 
distribution  with  zero-mean  and  covariance  ol.  (Eor  more  information  on 
circular  Gaussian  distributions  see  [20]).  Therefore,  random  matrix  Ec  which  is 
formed  from  N  random  vectors  and  given  by  Equation  (24),  is  circular  Gaussian 
distributed  (subscript  "c"  denotes  circular): 
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Ec  =[e(l),e(2),  -  ,e(N)  ; 


(24) 


Since  e(t)  =  e(t)  +  je(t)  and  both  e(t)  and  e(t)  are  real-valued,  a  real  random  matrix 

Er  can  be  formed  as  in  Equation  (25),  and  each  random  vector,  either  e(t)  or  e(t) 

forming  Er,  has  a  Gaussian  distribution  with  zero  mean  and  covariance  of  (o/2)I 
(subscript  "R"  denotes  real); 

=[E,E]  =  [e(l),  1(2),  •••  ,  e(N),e(l),e(2),  •••  ,  e(N)]  (25) 

Since  the  real  and  the  imaginary  parts  are  statistically  independent,  the 
multidimensional  Gaussian  pdf  [21]  for  the  measured  data  y(t)  is  dependent  on 
the  set  of  unknown  parameters  0  and  can  be  written  as  in  Equation  (26): 

F.(Y)  =  nP(X(t)  £  y(t),e)  =  nP(Y(t)+jY(t)  <  y(t)+jy(t),e) 

t-1  t-1 

=  nP(Y(t)  <  y(t),  e)nP(Y(t)  <  y(t),e)  =  F.(Y)F.(Y) 

t-1  t-1 

fe  (Y)  =  V,F,  (Y)  =  V,F,  (y)  x  V,F,  (y)  =  f,  (y)  x  (y) 

fe  (Y)  =  (27i)’""^  {det(K)}'^"^  exp|-^^(y(t)-Ax(t))^  K'*  (y(t)-Ax(t))| 
x(2ji)'""^{det(K)}'^exp|-^^(y(t)-Ax(t))^K-‘(y(t)-Ax(t))| 

Since  covariance  matrix  K  is  given  as  (o/2)I  for  each  y(t)  and  y(t).  Equation  (26) 
takes  the  form  of: 
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/  \-Nra  f  1  /  N  N  ^ 

fe  (Y)  =  (27i)'''”  ^  exp  —  ^(y(t)-Ax(t))^  (y(t)-Ax(t))  +  ^(y(t)-Ax(t))^  (y(t)-Ax(t)) 


/  \-Nm  f  1  N  I 

fjY)  =  (27if"[^|j  exp|--|j(y(t)-Ax(t)f  (y(t)-Ax(t))|  (27) 


The  likelihood  and  log-likelihood  functions  simply  become. 


/  \-Nm  r  ^ 

1(0,Y)  =  (27i)'"'"  ^  exp— |]e(t)“  e(t) 

V2;  [a  ,=i 


(28) 


1  N 

L(0,Y)=ln{l(0,Y)}=const -Nm  ln(a)  — (29) 


a  t=i 


where  the  vector  of  unknown  parameters  0  =  [o  x^(t)  x^(t)  0r]T 

The  derivation  of  CRB  as  follows  from  [6].  The  derivate  of  the  log- 
likelihood  function  L(0,Y)  is  taken  with  respect  to  the  unknown  parameters  o, 
x^(t),  x^(t),  0R  in  Equations  (30),  (31),  (32),  and  (33): 

s(a,  y)  =  |-  l(0,  y)  =  - —  +  4  Z  e“  (t )  e(t )  (30) 

oa  a  a  721^ 


s(x(k),  Y )  =  l(0,  y)  =  -^  [a“  e(k)+  A^  e*  (k)] 

r\ 

=  -Re{A»e(k)}  k  =  l,---,N 


(31) 


s(— ),Y)  =  ^L(0,Y)  =  i 
=— Im{A“e(k)} 


jA«e(k)+jA^e*(k)] 

k  =  l,---,N 


(32) 
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(33) 


s(0R.,Y)=^L(0,Y)  =  ^X*^4^“(t)^5(t)| 


ae 


Ri 


a 


t=i 


Ri 


i  =  l,---,n 


Define  new  matrices,  D  and  X  as: 


D  =  [d(0,J,-,d(0,J] 


da(0Ri)  da(0R^) 
d0R,  ’■■■’  d0R„ 


X(t)=diag(x(t))  = 


,(t)  0 


X 


„(t) 


Using  D  and  X,  Equation  (33)  can  be  rewritten  as: 


N 


s(e„  Y)  =  ^L(e,  Y)  =  -  X;  Re{x"  (t)D"  e(t)} 

d0R  a 


(34) 


CRB  AS  is  given  by: 

J-'={e[s(0,y)s"(0,y)]}‘‘  (35) 

where  s(0,Y)  is  called  the  score  function,  and  J  is  the  Fisher  Information  matrix  as 
explained  in  Section  2.3.  The  score  function  is  given  by  Equation  (10). 

The  derivation  of  the  Eisher  Information  matrix  is  shown  below.  The 
covariance  of  s(^Y)  must  be  computed.  The  first  step  is  to  compute  the 
correlation  of  each  unknown  parameter  with  itself,  and  with  the  other  unknown 
parameters. 
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i(a,Yf 


mN  1  H  /■  \  / 
- +  ^^5  (t)e(t) 


a  a  73 


t=i 


) 


2mN  ■  ^ 


N  N 
^  t=l  s=l 


X  (t)  e(t)]+  ^  S  Z  (t)  e(t)  e“  (s)  e(s) 


Since  E 


(t)  e(t)  =ma  and  e(t)  independent  of  e(s)  for  all  s/t 


e(t)  e“(s)  e(s) 


2  2 
im  a 


for  t  s 


m 


(m+l)a^  fort  =  s 


(for  proof  see  PRF  1  in  Appendix  A),  E[s(o,Y)2]  takes  the  form  of: 


i(a,Y)" 


2m"N"  Nm 


a 


2  2 

a  a 


+  ^[(N-l)m+(m  +  l)]  = 


mN 


a 


(36) 


Using  E  e”(t)e(t)e^(s)=0  (See  Appendix  A  for  proof),  s(o,Y)  is 

uncorrelated  with  the  other  score  functions  for  the  remaining  unknown 
parameters. 

To  simplify  the  equations  written  below,  transformations  between  the  real 
and  imaginary  parts  of  different  complex  variables,  which  are  given  in  Appendix 
A  for  all  t  and  s,  are  used. 


E  s(x(k),Y)  (x(p),Y)  is  found  below  using  (31). 


E  s(x(k),Y)  (x(p),Y)  =E  — Re|A“e(k)|  — Reje^  (p) A*| 

L  J  r7  a  ^ 


4  1 


2  ^|Re|A“E  e(k)  e^(p)  A*|+Re|A^E  e(k)e^(p)  Ajj 


(37) 


e(t)  e^(s) 


=  0 


(38) 


40 


E  e(t)  (s)  =  al 


(39) 


Substituting  Equations  (38)  and(39)  in  Equation  (37)  leads  to  Equation  (40). 

E[s(x(k),Y)s'(l(p),Y)]  =  -Re{A"A}8^,  (40) 

a 

E[^s(x(k),  Y)s^  (x(p)?Y)J  is  found  in  a  similar  manner  from  Equations  (31),  and 
(32). 

"(p)A*} 

^  (41) 

=  ^(lm{A“E  e(k)e^(p)]A*}-Im{A“E  e(k)e“(p)]A} ) 
a 

Substituting  Equations  (38)  and(39)  in  Equation  (41)  leads  to  Equation  (42). 

E[s(x(k),Y)sMx(p),Y)]=  --Im{A"A}8t,  (42) 

a 

E[^s(x(k),Y)s^  (Qr^Y)  is  found  below  using  Equations  (31),  and  (34). 

E[s(x(k),  y)s"  (0J, ,  y)]  =  Er4Re{A«  e(k)}X  Reje"  (t)D*X*  (t)} 

""  (43) 

=  4  Z  (Re{A«E[e(k)e"  (t)]D‘X*  (t)}+  Re{A«E[e(k)e“  (t)]DX(t)}) 

a  t-i 

Once  again  using  Equations  (38)  and  (39),  Equation  (43)  simplifies  to  Equation 
(44). 

E[s(x(k),  y)s"  (0  ,  y)]  =  4  Re{A«DX(k)}  (44) 

a 

E[^s(x(k),Y)  (x(k),Y)J,  and  E[^s(x(k),Y)  (0r,Y)J,  which  are 

obtained  in  the  same  manner  that  E[^s(x(k),Y)  (x(k),Y)J,  and 


E[s(x(k),Y)s^(x(p),Y)]=  E  — Re{A“e(k)}— Imje 

a  a 
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E  s(x(k),Y)  (0r,Y)  are  obtained  in  Equations  (40),  and  (44),  are  given  below 


in  Equations  (45),  and  (46). 


E[s(x(k),  Y)sMx(p),  Y)]=  4Re{A»E[#)e"(p)]A}=  (45) 

a  a 


E[s(x(k),Y)s^(0R,Y)]  =  --Im{x“(k)D“Af  = -ImjA^^DXlk)}  (46) 


E  s  (0R  ,Y )  (0R  ,Y )  is  shown  as  the  last  step. 


e[s(0^,y)s^(0„y)]  =  E  4ZZRe{x«(t)D«e(t)}Reie"(s)D*X*(s)} 


4 

O'  t-1  s-1 


(Re{x"  (t)D“E[e(t)e^  (s)]d*X*  (s)}+  Re{x"  (t)D“E[e(t)e“  (s)]dx(s)}) 


2 

O  t=l  s=l 


Substituting  Equations  (38)  and  (39)  leads  to  Equation  (47). 

E[s(e,,Y)s''(e,,Y)]=-|;Re{x“(t)D“DX(t)} 


(47) 


Now  the  complete  Fisher  Information  matrix  can  be  constructed: 

e[4».y)"]  0 

£[qr(l).Y)  /  b(l).Y)]  i;[4x(l).Y)  /  (i(l).Y)] 

£[,(1(1). y)  7 (x(i).y)]  £[440. y) 7 (i(0.d] 


£[4^(1), Y)  7  (0.,Y)] 
£[4y(1),Y)  7(0„Y)] 


(£[.,(r(l).Y). y)])  (£[.,(i(l).Y). y)]) 


£[4x(y).y)  7(x(jv).y)]  £[4x(q.Y)  7(i(w).Y)]  £[£(x(V),Y)  7(^„,Y)] 
£[4i(iv).Y)  7  (x(w).y)]  £[.4i(iv).Y)  7  (i(jv).Y)]  £[£(x(V),Y)  i''  (^j.Y)] 

(£[.,(j(JV).Y)7(£„,y)])  (£[.,(i(Y).Y)7(s„,Y)])  £[£(e,,Y)  /((S^.Y)] 


(48) 


The  following  terms  are  introduced  for  notational  convenience. 
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-yRe{x“(t)D“DX(t)} 

a  t=i 


var. 


CR 


(a)  = 


a 

mN 


H  =  -A“A 
a 

G  =  H  ‘ 

Ak=-A“DX(k) 

a 


G  is  defined  with  respect  to  H  and  they  are  related  in  an  inverse  sense. 
This  special  relation  between  H  and  G  becomes  useful  while  inverting  the  Fisher 

A 

Information  matrix.  Simply  G=H“',  but  the  following  relation  also  holds  (for 
proof  see  PRF  3  in  Appendix  A). 


H 

-h" 

-1 

G 

-g" 

H 

hJ 

G 

G 

The  H  matrix  is  Hermitian  and  therefore  its  imaginary  part  must  be  skew- 
symmetric  =  -H . 

J(0)  introduced  in  Equation  (48)  can  be  rewritten  by  using  the  newly 
defined  terms: 


var-^(a) 


H  -H 
H  H 


0 


H 

-H 

An 

H 

H 

An 

XT  XT 

All  ^1 

•• 

r 

(49) 
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In  order  to  obtain  the  CRB  as  a  function  of  DOA  angles  (0r),  J(0r)  is 
needed.  J(0r)  is  found  in  (50)  by  partitioning  J'^. 


j(0j  =  r-p Af  ...a;; a^] 


G  -G 

"a  " 

0 

1 

G  G 

G  -G 

0 

G  G 

A_ 

(50) 


Using 


G  -G 
G  G 


G  A-GA 
G  A  +  GA 


Re(GA) 

Im(GA) 


and 


Re(GA) 

Im(GA) 


Equation  (50)  takes  the  form  in  Equation  (51). 


j(0R)-r-|;Re[A»GA] 

2  f  .  ^  .  1 

J(0  R )  =  -  Z  Re  [x  ”  (t)D  “  DX(t)  -  X  “  (t )d  “  a(a  ”  a)  A  “  DX(t) 

Finally  in  Equation  (52)  J(0r)  is  inverted  in  order  to  obtain  the  formula  for 
CRBas(0r). 


CRB^,(0r)-J'‘(0r) 


f  N 

ZRe 

X^(t)D“ 

I-A(a“a)‘‘ A^ 

DX(t) 

Y 

V  t=i 

- 

- 

- 

/ 

(52) 


For  sufficiently  large  N,  the  CRBas(0r)  is  given  by  Equation  (53)  (see 


Appendix  B  for  proof).  The  dependence  of  the  CRB  as  on  0r  is  shown: 

CRB„  (e  J  =  [r4(d-  [l  -  a(a"  a)-'  A"  ]d).  *?'}['  (53) 


2N 


where  P  =  E 


x(t)  x^(t)  Aod  "  is  array  multiplication  operator. 
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Equation  (53)  implies  that  when  N  is  sufficiently  large,  the  CRB  can  be 


determined  depending  on  o,  P,  and  the  true  DOA  of  the  signal  i  (0Ri)  without  the 
need  for  the  eigenvalue  decomposition  of  the  covariance  matrix  R  of  the 
observation  data,  y(t).  In  this  case,  the  CRB  for  0Ri  is  given  by  (54). 

CRB„(eK.)=^W(D-[l-A(A"ArA"]D).*P'}]^'  (54) 


In  our  model,  the  number  of  realizations  N  is  512  which  makes  Equation 
(54)  a  good  approximation  for  determining  the  CRB  as  for  different  DOA  and 
different  number  of  signals  with  different  correlation  levels. 


3.3  Derivation  of  the  CRB  for  Constant  Modulus  Signals 
(CRBcm) 

Constant  Modulus  signals  have  constant  envelope  and  these  signals  are 
modulated  via  a  change  in  frequency.  Some  examples  of  constant  modulus 
signals  include  Erequency  Modulation  (EM),  Phase  Modulation  (PM),  Erequency 
Shift  Keying  (ESK),  and  Phase  Shift  Keying  (PSK).  The  Constant  Modulus 
Algorithm  (CM A)  was  first  considered  in  [22,  23]. 

Derivation  of  CRBcm  follows  [24].  Eor  CM  signals,  measured  data  is 
modeled  by, 

y(t)=ABs(t)+e(t)  (55) 

where 

y(t)  is  m  by  1  measured  data  vector  at  time  t. 
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A  =  [a(0Ri),---,a(0j^n)]  where  a(0Ri)  is  m-by-1  array  response  vector 
from  direction  0Ri,  and  0r  =  [0ri,  . . . ,  0Rn]  is  the  DOA  vector  of  the 
signals, 

B  is  n-by-n  channel  gain  matrix  and  has  the  amplitudes  of  the 
signals  =  [|3i,  ,  (3n])  on  the  main  diagonal.  Amplitude  of  the 

i*  signal  as  received  by  the  antenna  is  a  real  parameter  and  is 
denoted  as  |3i, 

s(t)  is  an  n-by-1  vector  of  source  signals  at  time  t  with  unit 
amplitude,  equal  to[si(t),  ... ,  Sn(t)], 
e(t)  is  an  m-by-1  additive  noise  vector  with  complex  parameters. 

Different  signal  amplitudes  are  absorbed  into  the  B  matrix.  Therefore  the 
s(t)  vector  only  has  phase  information  corresponding  to  observation  time  t. 
Hence,  the  phase  of  the  i*  signal  can  be  written  as  S;  (t)  =  where  O;  (t)  is  the 
unknown  phase  modulation  for  i*  signal.  A  source  vector  can  be  defined  as 
0(t)  =  [cl)j  (t),  •••  (t)]^,  the  phase  vector  for  all  signals  at  time  t. 

Maintaining  assumptions  of  additive  noise  as  for  the  AS  case,  the 
likelihood  function  given  in  Equation  (28)  also  holds  true  for  CM  signals  except 
that  the  form  of  e(t)  is  different.  In  each  case  Likelihood  and  log-likelihood 
functions  are  given  in  Equations  (56),  and  (57)  respectively: 

/  \-Nm  f  1  N  t 

l(0,Y)  =  (27i)'''"'f|j  exp|--|](  y(t)“-s“(t)  B'^A")  (  y(t)  -ABs(t)  (56) 
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L(0,Y)=const-Nmln(a)-  — ^fy(t)”  -  (t)  j  (y (t)  -  ABs(t))  (57) 

a  ,=i  “ 


Taking  derivates  of  the  log-likelihood  function  with  respect  to  the 
unknown  parameters  ®(t),  ^(t),  and  0r  leads  to  the  individual  score  functions. 


The  set  of  unknown  parameters  are  represented  by  a  vector,  denoted  as  ^  where 


0= 


0"(1),  •••  ,0"(N),P"(1),  •••  ,P"(N),0 


The  partial  derivative  with  respect  to  0(t)is  found  by. 


and  _  'L(0,Y),  ^  L(9,Y)  follows  from  Equations  (31),  and  (32). 

3s(t)  os(t) 


^L(e,Y)=iRe{B"A"e(t)} 
j^L(e,Y)  =  hm{B"A"s(t)} 


Since  the  s(t)  vector  only  carries  phase  information  (the  observed  signal 
amplitudes  are  contained  in  the  matrix  B),  its  real  and  the  imaginary  parts  are 


s(t)  =  Re{s(t)}  =  cos($(t)),  and  ^(t)  =  Im-^(t)}=  sin($(t)).  Therefore 


^s(t) 

3$(t) 


and 


^s(t) 

3$(t) 


are  simply  found  as: 


30(t) 


=  -sin(o(t)) 


’  5®(t) 


=  cos(o(t)) 
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Let  Xt  be  a  matrix  that  has  s(t)  on  its  main  diagonal,  =diag(s(t));  then 


^s(t)  ,  dljt) 

a#)'  a#) 


can  be  written  in  terms  of  St. 


^s(t) 

a$(t) 

^s(t) 

a$(t) 


Im{sJ 


Re{sJ 


Then  Equation  (58)  becomes, 

s($(k),Y)  =  ^^L(0,Y)  =  ^(-Im{xjRe{B«A“e(k)}+Re{xJlm{B«A»e(k)}) 


s($(t),Y)  =  -Im{sfB“A“e(t)} 
a 

s(0R ,  y)  is  found  with  a  similar  manner  as  Equation  (34) 


s(e„Y) 


=  ^L(e,Y)=-j;Re{s:<B"D"e(t)} 

a0R  a 


where  D  =  [d(0Ri),---,d(0R„)] 


da(0Ri) 

d0Ri 


s((3,y)  is  obtained  from  [6]. 


s(p,  y)=  4  L(e,  Y)  =  - 1;  Re{s>"  e(t)} 

-  ap  a 

Fisher  Information  matrix  can  be  derived  as. 


(59) 


(60) 


(61) 
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Hi  0  E;' 

J=  0  El  (62) 

Ai  •••  A^  r 

Ej  •  •  •  Ej^  A  M 

where 

Hi=E[s($(t),Y)s^($(t),Y)]=-Re{sfB»A»ABSi}  (63) 

a 

Ai  =e[s(0i,,y)s^(#),y)]=  --Im{sfB«D“ABSi}  (64) 

a 

Ei=E[s(p,Y)s^(#),Y)]  =  --Im{s»A«ABSi}  (65) 

a 

r  =  E[s(0i,,Y)s"(0i,,Y)]=--XRe{srB«D«DBS,}  (66) 

a 

A=E[s(p,Y)s"(0i,,Y)]=--XRe{s«A«DBSi}  (67) 

a  t::rl 

M=E[s(p,Y)s"(p,Y)]=--|;Re{sfA«ASi}  (68) 

<5  t=i 


Inversion  of  the  Fisher  Information  matrix  is  carried  out  in  block- 
partitioned  form  using  Schur  complement  formulas  and  the  Woodbury  identity 
as  in  [24],  The  CRB  for  DOA  of  CM  signals  is  given  by: 

CRB„(eJ=diag(>P„  (69) 

where 


Til  A^ 

'*'21  '*'22]  am 
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B„=£a,HX 

t^l 

Bj,  =Ze,h;'a; 

t-1 

b„=Ea,hX 

t  =  l 

B2=i;E,H;'E; 


t=l 


3.4  CRB  Calculation 

The  CRBs  for  different  scenarios  are  calculated  and  plotted  by  using 
Matlab®  software  to  illustrate  performance  under  different  conditions.  The  first 
scenario  is  where  only  one  object  is  assumed  to  exist  in  the  surveillance  region. 
The  CRB  for  different  object  locations  in  azimuth  (0r)  and  different  SNR  levels  is 
observed. 

The  second  scenario  includes  two  objects.  One  of  the  objects  simulates  a 
moving  object  and  its  DOA  changes  between  Coherent  Processing  Intervals 
(CPIs)  and  as  it  moves  from  one  angular  location  in  azimuth  to  another.  The 
second  object  remains  stationary  at  a  fixed  location  in  azimuth  throughout  all 
CPIs.  A  stationary  object  does  not  necessarily  imply  the  object  is  a  non¬ 
maneuvering  object;  it  implies  the  object  does  not  change  location  in  azimuth. 
The  CRB  for  both  objects  is  plotted  separately  on  different  figures. 

Interference  from  multiple  transmitters  in  the  surveillance  region  is  one  of 
issue  in  PCL  systems.  Thus,  to  observe  the  effects  of  having  a  transmitter  in  the 
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surveillance  region  on  the  bounds  of  the  object  echo  signals,  a  third  scenario  is 
introduced  and  includes  two  object  echo  signals  and  a  direct  signal  from  a 
transmitter  of  opportunity  in  the  surveillance  area.  The  transmitter  is  simulated 
by  choosing  a  relatively  high  SNR  level  compared  to  the  SNR  levels  of  the  object 
signals,  and  having  zero  Doppler  shift.  The  object  signals  are  chosen  to  be  a 
moving  object  and  a  stationary  object  as  in  the  second  scenario.  Last  but  not 
least,  fourth  scenario  observes  the  CRB  of  the  moving  target  for  AS  and  CM  cases 
for  different  transmitter  of  opportunity  DOAs.  This  experiment  should  yield 
additional  insight  into  antenna  placement  and  orientation. 

Correlation  between  object  returns  is  also  implemented  for  the  arbitrary 
signal  case.  The  matrix  P,  defined  in  Equation  (70),  includes  the  correlation  of 
the  signals.  In  the  case  when  two  objects  have  the  same  SNR  values  and  they  are 
uncorrelated,  the  correlation  coefficient,  p-0.  For  fully  correlated  signals,  p  -1. 


P  = 


1 

P 


P 

1 


(70) 


SNR  variation  is  implemented  by  assigning  new  values  of  o.  Signal 
power  is  assumed  to  be  1  Watt  and  the  relation  between  the  o  and  the  SNR  level 
is  given  in  Equation  (71)  for  SNR  expressed  in  dB: 


a  = 


1 


(71) 
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Dimensions  of  P  change  according  to  the  number  of  signals  being 
received.  To  illustrate  how  the  cross  correlation  of  multiple  signals  for  different 
numeric  SNR  levels  take  place,  a  matrix  P  for  three  signals  is  given  in  Equation 
(72)  (See  Appendix  C  for  proof).  Signal  one  has  been  normalized  to  1  watt. 


3.5  Summary 

This  chapter  explained  and  showed  the  derivation  of  CRB  as  and  CRBcm. 
The  implementation  of  correlation  between  two  signals  is  also  explained  and 
expanded  to  multiple  correlated  signals  for  AS  case.  The  formulas  given  here  are 
used  for  the  CRB  calculations. 

Equation  (54)  which  is  approximated  for  sufficiently  large  number  of 
realizations  (N)  is  used  for  the  CRB  as  calculations. 

CRB„(e„)  =  ^[Re{(D-[l-A(A"ArA"]D).*P"}]:'  (54) 

Equation  (69)  is  used  for  the  CRBcm  calculations. 

CRB«(e,)=diag(>F„  (69) 
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4  .  CONTRIBUTIONS  OF  RESEARCH 


This  chapter  presents  results  computed  for  the  CRB  of  Arbitrary  Signals 
and  Constant  Modulus  signals  in  different  scenarios.  The  CRBs  for  moving 
signal  sources  are  plotted  versus  the  true  DOA  of  the  moving  signal  sources,  and 
CRBs  of  the  stationary  signal  sources  are  plotted  versus  angular  separation  of  the 
stationary  signal  sources  with  the  moving  signal  sources.  Matlab®  software  is 
used  for  calculating  the  CRBs  and  for  generating  the  plots. 

4.1  Scenario  1:  Single  Signal 

The  simplest  case  considered  includes  only  one  signal  impinging  on  the 
antenna.  This  case  is  studied  because  it  is  a  good  way  of  evaluating  the  system 
performance  before  moving  on  to  more  complicated  scenarios.  The  single  signal 
provides  an  understanding  of  basic  system  performance  and  gives  a  reference  for 
comparison  when  studying  more  complicated  scenarios. 

CRBas  and  CRBcm  are  equal  to  each  other  when  there  is  only  one  signal. 
This  is  because  the  CM  signal  model  in  Equation  (55)  can  also  be  represented 
using  the  arbitrary  signal  model  given  in  Equation  (21)  in  the  single  signal  case. 
Since  the  pdfs  of  CM  and  arbitrary  signals  are  basically  identical  (differing  only 
in  the  way  that  their  means  are  represented),  then  the  single  signal  case  results  in 
identical  pdfs  for  the  CM  and  arbitrary  signals. 
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Figure  12  shows  the  CRB  of  a  single  object.  The  CRB  is  a  function  of  the 
true  DOA  of  the  signal,  or  in  other  words,  it  is  a  function  of  the  angular 
difference  between  object  location  in  azimuth  and  the  antenna  broadside 
direction.  As  the  location  of  the  signal  source  in  azimuth  shifts  away  from  the 
antenna  broadside  direction,  there  is  a  nonlinear  increase  in  the  CRB.  The  CRB  is 
symmetric  in  this  case  because  it  is  a  function  of  the  signal  angular  location  in 
azimuth,  as  measured  from  the  antenna  broadside  direction.  For  example,  a 
signal  source  has  the  same  CRB  when  it  is  either  located  at  -10°  or  10°  from 
antenna  broadside  (assuming  all  other  parameters  such  as  SNR  remain  constant). 


SNR  =  0dB 


Figure  12.  CRB  for  single  signal  case  when  SNR  =  OdB. 

An  increase  or  decrease  in  SNR  results  in  a  linear  change  on  the  CRB  in  an 
inverse  sense;  increasing  SNR  causes  a  decrease  in  the  CRB  and  decreasing  SNR 
causes  an  increase  in  the  CRB.  The  effect  of  SNR  changes  on  the  CRB  is 
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demonstrated  in  Appendix  D  where  the  CRBs  for  -lOdB,  lOdB,  and  20dB  SNRs 
are  plotted.  The  shape  of  the  CRB  does  not  change  but  there  is  an  increase  or  a 
decrease  of  the  CRB  amplitude  inversely  proportional  to  the  change  in  the  SNR. 
This  result  is  expected;  as  signal  power  increases  one  expects  the  variance  in  the 
DOA  estimate  to  decrease. 

4.2  Scenario  2:  Two  Signal  Sources 

When  receiving  multiple  one  signals,  the  CRBas  and  CRBcm  show 
differing  characteristics.  This  is  due  to  the  CM  constraint  on  the  received  signal 
type  [25]. 

DOA  estimators  that  utilize  signals  with  the  CM  property  separate  the 
incoming  signals  using  different  techniques,  and  then  utilize  the  signal 
properties,  such  as  signal  phase,  in  addition  to  the  array  antenna  manifold,  that 
conventional  DOA  estimators  use.  Not  only  is  the  phase  difference  between 
array  antenna  elements  used  to  estimate  DOA,  but  also  the  actual  phase 
information  of  the  incoming  signal. 

Since  the  phase  of  the  incoming  signal  is  also  processed  in  CM  estimators, 
a  derivative  is  taken  with  respect  to  the  signal  phase  when  the  Fisher  Information 
matrix  is  being  constructed.  The  signal  phase  is  added  to  the  vector  of  unknown 
parameters  when  the  CM  constraint  is  brought  into  the  estimator.  This  is  the 
reason  why  CM  estimators  have  a  different  CRB  than  AS  estimators. 
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4.2.1  CRB  AS  of  Two  Signal  Sources 

In  Figure  13,  the  CRB  for  the  so-called  "moving  object"  is  illustrated.  One 
signal  is  a  stationary  signal  source  sitting  at  0°  and  the  other  signal  source  is  a 
moving  signal  source.  It  is  called  a  "moving  signal"  because  the  source  remains 
fixed  at  one  angle  during  512  realizations  of  measurement  data,  after  which  it  is 
moved  to  the  next  DOA  where  it  remains  fixed  for  another  512  snapshots. 
Basically,  the  moving  object  remains  fixed  during  a  CPI  but  is  moved  to  the  next 
DOA  angle  between  CPIs.  The  stationary  signal  source  remains  at  the  same  exact 
location  throughout  all  CPIs.  This  is  a  fair  approximation  to  a  moving  object 
scenario;  recall  from  Section  3.2  that  the  number  of  realizations  must  be  kept 
large  (in  this  case,  N  =  512)  for  the  CRB  approximation  of  Equation  (53)  to  remain 
valid. 

When  CRB  AS  plots  are  observed  it  can  be  seen  that  the  CRB  as  out  to  70°  is 
relatively  high  due  to  the  system  performance  shown  in  single  signal  case  (see 
Figure  13  (b)).  This  is  actually  seen  in  the  CRB  plots  for  single  signal  case.  As  the 
moving  object  gets  closer  to  the  other  signal  source,  a  nonlinear  increase  in  the 
CRB  AS  is  observed  (see  Figure  13  (a)  and  (c)).  This  is  caused  by  both  signals 
having  close  DOA  angles,  thus  causing  the  phase  delays  between  the  elements  of 
the  array  antenna  to  become  very  similar.  Hence,  the  array  response  vectors 
a(0i),  and  a(02)  become  more  alike  as  the  angular  difference  between  the  two 
signals  is  reduced  in  the  azimuth  plane.  When  two  signals  are  exactly  on  the  top 
of  each  other,  the  A  matrix  formed  by  a(0i)  and  a(02)  and  the  D  matrix  formed  by 
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-^a(0i)  and  Tr^a(02)A  become  singular  and  cause  the  CRB  as  to  approach 

C/O  j  C/O  2 


infinity.  Therefore,  as  the  DOA  of  two  signals  gets  close  to  each  other  the  CRB  as 
increases.  When  they  are  exactly  at  the  same  DOA  we  say  that  it  is  unbounded. 


Moving  Target  SNR  =  OdB  Moving  Target  SNR  =  OdB 


(a)  (b) 

Moving  Target  SNR  =  OdB 


(c) 


Figure  13.  The  CRB  as  of  a  moving  object  with  OdB  SNR  with  a  stationary  object  of 

the  same  SNR  located  at  0°. 
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Figure  13  (b)  shows  the  decrease  in  CRBas  as  the  moving  signal  source 
moves  away  from  70°  off  broadside  towards  broadside.  Figure  13  (a)  shows  the 
increase  in  CRBas  in  the  DOA  angle  range  of  1°-10°.  When  the  moving  object  is 
at  1°  of  DOA  its  CRBas  is  nearly  ten  times  as  more  than  when  it  is  out  at  70°. 
Finally,  the  DOA  range  of  0°  -  1°  is  shown  in  Figure  13  (c)  and  the  increase  in 
CRBas  is  shown  when  two  signal  sources  are  very  near  to  each  other. 

To  gain  a  better  understanding  of  the  increase  in  CRBas  as  the  angular 
separation  between  the  signals  decrease,  the  CRBas  of  the  stationary  signal 
source  is  plotted  versus  its  angular  separation  with  the  moving  object  in  Figure 
14.  Therefore,  the  CRBas  of  the  stationary  object  is  calculated  for  each  DOA  of 
the  moving  signal  and  plotted  versus  the  angular  separation  with  the  moving 
signal  source. 

The  CRBas  of  the  stationary  object  does  not  exhibit  a  monotonic  increase 
or  decrease.  In  Figure  14  (a),  the  minimum  CRBas  of  the  stationary  object  is 
achieved  when  the  moving  signal  source  is  furthest  away.  As  the  two  signal 
sources  get  closer  in  azimuth,  the  difference  between  the  minimum  and 
maximum  points  increases.  The  CRBas  values  in  Figure  14  (b)  and  (c)  are  similar 
to  CRBas  values  of  the  moving  object;  this  is  because  the  A  and  D  matrices 
become  nearly  singular  as  two  signal  sources  gets  closer  to  each  other  in 
azimuth. 

Correlation  between  the  arbitrary  signals  should  also  be  incorporated  into 
our  DOA  estimate,  since  the  bounds  change  with  respect  to  the  correlation  level 
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of  the  signals.  The  effect  of  correlation  between  two  signals  is  shown  in 
Appendix  E,  where  bounds  for  the  moving  and  the  stationary  signal  sources  are 
plotted.  The  higher  the  correlation  between  the  arbitrary  signals,  the  higher  the 
bound  becomes. 


(a) 


(b) 


stationary  Target  SNR  =  OdB 


(c) 


Figure  14.  The  CRB  as  of  a  stationary  object  with  OdB  SNR  with  a  moving  object 


of  the  same  SNR. 
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When  the  CRB  as  in  the  two  signals  case  is  compared  with  the  CRB  in 
single  signal  case,  the  bound  increases  in  the  antenna  broadside  direction  in  the 
CRB  AS  of  two  signals  when  the  DOA  of  the  two  signals  are  near  to  each  other. 
4.2.2  CRBcm  of  Two  Signal  Sources 

As  mentioned  in  Section  3.3,  the  CM  constraint  on  the  signal  changes  the 
CRB.  In  this  research  the  CRBcm  is  unique  because  of  the  effects  of  the  changes 
in  the  Doppler  shift  and  the  start  phase  of  the  signal  sources.  These  parameters 
are  modeled  using  the  FM  signal  model.  Estimators  having  the  CM  signal 
constraint  use  many  signal  properties,  especially  properties  of  the  signal  phase. 
Therefore  different  effects  on  the  phase  of  the  incoming  signals,  such  as  different 
Doppler,  different  start  phases  of  the  signals,  or  different  information  messages 
being  modulated  can  cause  the  signals  to  be  distinguished  from  one  another. 

As  different  signal  sources  get  closer  to  each  other  in  azimuth,  the  A  and 
D  matrices  become  singular,  but  the  CRBcm  is  not  as  susceptible  to  reduced 
angular  separation  between  signal  sources  in  azimuth  plane  as  the  CRB  as.  This 
is  because  of  the  CM  constraint  on  the  signal;  a  new  unknown  parameter  vector 
of  phase  must  be  added  to  the  unknown  parameter  vector  (0  for  arbitrary 
signals.  Hence,  the  derivative  of  log-likelihood  function  with  respect  to  phase 
vector  of  the  signal  sources  is  found  in  Equation  (59).  The  correlation  of  the 
unknown  parameter  vector  phase  with  itself  and  the  other  unknown  parameters 
given  in  Equations  (63),  (64),  and  (65)  makes  the  Eisher  information  matrix 
invertible  even  when  two  signal  sources  are  closely  located  in  the  azimuth  plane. 
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(C)  (d) 

Figure  15.  The  CRBcm  of  the  moving  signal  source  with  OdB  SNR  with  a 
stationary  signal  source  of  the  same  SNR  when  the  Start  Phases  and  the  Doppler 
Shifts  of  the  two  signal  sources  are  the  same. 

The  CRBcm  is  also  found  to  be  a  function  of  start  phase  and  Doppler  shift 
differences.  In  Figure  15,  the  start  phase  and  the  Doppler  shift  of  the  incoming 
signals  are  the  same,  so  this  can  be  thought  of  the  case  when  the  phases  of  the 
two  incoming  signals  are  fully  correlated.  Due  to  the  reason,  the  signal  model 
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used  in  this  thesis  assumes  identical  FM  modulation  (message)  signals  for  the 
two  signals,  and  that  this  FM  signal  remains  constant  during  the  sampling  of  the 
data.  Therefore,  having  the  same  Doppler  shift  and  starting  phase  for  the  two 
signal  sources  should  be  considered  as  a  worst-case  scenario. 

In  Figure  15  (b)  the  CRBcm  for  the  moving  object  seems  to  be  almost 
identical  to  the  CRBas  of  the  moving  object  within  the  DOA  range  of  10°  -  70°. 
The  real  problem  is  in  the  two  signals  case,  especially  when  the  two  signals  get 
closer  in  azimuth;  this  research  now  focuses  on  the  region  where  the  signals  are 
relatively  close  in  azimuth.  When  the  moving  object  is  within  10°  of  the 
stationary  object,  its  CRBcm  is  initially  lower  than  CRBas.  But,  as  the  moving 
object  approaches  the  stationary  object  at  0°,  the  increase  in  the  CRBcm  is  greater 
than  the  increase  in  CRBas.  This  is  because  to  the  two  signals  have  exactly  the 
same  phase  properties  (due  to  our  initial  assumptions  of  identical  starting  phase 
and  Doppler). 

When  the  two  signals  have  a  starting  phase  difference  and  a  Doppler  shift 
difference,  the  CRBcm  varies  because  it  can  better  separate  two  signals  from  each 
other.  This  is  shown  in  Figure  16  by  choosing  an  arbitrary  starting  phase  and 
Doppler  shift  difference  between  the  signals. 

Figure  16  (a)  shows  that  when  there  is  a  difference  in  the  phase 
parameters  of  the  two  signal  sources,  the  CRBcm  of  the  moving  object  remains 
the  same  as  the  CRBas  of  the  moving  object  in  Figure  13  (a)  and  CRBcm  of  the 
moving  object  in  Figure  15  (a)  between  the  DOA  range  of  10°  -  70°.  The  CRBcm 
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(a)  (b) 


(c) 


Figure  16.  The  CRBcm  of  moving  signal  source  with  OdB  SNR  with  a  stationary 
signal  source  of  the  same  SNR  when  Start  Phases/ Doppler  Shifts  of  the  moving, 
and  stationary  signal  sources  are  [2,4]/ [8,-2]  respectively, 
in  Figure  16  (a)  takes  a  better  look  at  the  behavior  of  the  bound  for  the  moving 
object  when  the  two  signals  approach  each  other  in  azimuth  plane.  CRBcm  does 
not  increase  as  much  when  compared  with  the  previous  cases  and  it  is  the  most 
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robust  bound  obtained,  considering  the  increase  in  the  CRB  of  the  previous  cases 
(when  the  two  signals  became  near  to  each  other  in  azimuth). 


Angular  Separation  (in  degrees)  Angular  Separation  (in  degrees) 


(a) 


(b) 


(c) 


(d) 


Figure  17.  The  CRBcm  of  the  stationary  signal  source  with  OdB  SNR  with  a 
moving  signal  source  of  the  same  SNR  when  the  Start  Phases,  and  the  Doppler 
Shifts  of  the  two  signal  sources  are  the  same. 
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The  final  CRBcm  analysis  for  this  case  involves  the  CRBcm  of  the 
stationary  signal  source.  Figure  17  shows  the  worst  case  scenario  for  the  CM 
signal  estimators  when  the  phase  parameters  of  the  two  signals  are  fully 
correlated.  As  a  result,  (just  as  in  the  case  of  the  CRBas  of  the  stationary  object) 
there  is  not  a  constant  increase  in  the  bound  as  the  separation  between  the  two 
signal  sources  gets  smaller.  After  noting  the  similarity  of  Figure  15  (a,  c,  d)  with 
Figure  17  (a,  c,  d),  it  becomes  apparent  that  once  the  separation  between  the 
moving  and  stationary  signal  sources  becomes  less  than  ten  degrees,  the  CRBcm 
of  the  stationary  object  and  moving  object  is  identical. 

As  a  result,  the  CRBcm  of  both  signal  sources  having  identical  phase 
parameters  yield  similar  bounds  with  the  CRBas  when  the  signal  separation  is 
between  10°  -  70°.  On  the  other  hand,  when  the  angular  separation  of  the  two 
signals  in  azimuth  is  less  than  ten  degrees,  the  CRBas  is  lower  than  the  CRBcm 
for  the  same  phase  parameters.  It  must  be  remembered  that  when  both  of  the 
signals  have  the  same  phase  parameters,  this  implies  both  signals  are  fully 
correlated  in  phase  and  represents  a  worst  case  scenario  for  CM  signals.  This  is 
rarely  the  case  in  real-world  scenarios. 

The  arbitrary  phase  parameters  used  for  generating  Figure  16  are  also 
used  to  generate  Figure  18.  When  the  separation  between  the  two  signal  sources 
are  in  the  10°  -  70°  range,  the  CRBas  behavior  illustrated  in  Figure  14  (a)  is  similar 
to  the  CRBcm  response  of  Figure  18  (b).  However,  when  the  signal  separation  is 
less  than  ten  degrees,  the  CRBcm  of  the  stationary  object  achieves  the  CRBcm  of 
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the  moving  object.  Note  that  the  plots  in  Figure  16  (b,  c)  for  the  moving  signal 
source  are  exactly  the  same  as  those  in  Figure  18  (b,  c)  for  the  stationary  signal 


source. 


(a) 


(b) 


(c) 


Figure  18.  The  CRBcm  of  the  stationary  signal  source  with  OdB  SNR  with  a 
moving  signal  source  of  the  same  SNR  when  Start  Phases/ Doppler  Shifts  of  the 


moving,  and  stationary  signal  sources  are  [2,4]/ [8,-2]  respectively. 
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The  CRBcm  of  signals  with  different  phase  parameters  seems  to  be  a  lower 
bound,  especially  when  angular  signal  separation  in  azimuth  gets  smaller.  When 
the  angular  separation  of  the  signals  in  azimuth  is  relatively  large  (in  the  10°  -  70° 
range),  the  CRBcm  of  the  signals  with  different  phase  parameters  seems  to  be 
slightly  lower  than  if  not  almost  the  same  as  the  CRB  as  of  two  signals  with  zero 
correlation. 

4.3  Scenario  3:  Two  Signal  Sources  and  a  Transmitter  of 
Opportunity  (Interferer) 

PCL  systems  can  be  susceptible  to  interfering  signals  in  the  environment, 
or  even  interfering  signals  from  their  own  transmitters  of  opportunity.  In  this 
section,  the  effects  of  a  transmitter  on  the  bounds  of  CM  and  arbitrary  signals  are 
illustrated.  Figure  19  shows  the  CRBas  when  all  three  signals  are  uncorrelated. 
Figure  20  shows  CRBcm  for  the  phase  parameters  that  are  specified  on  the 
figures. 

4.3.1  CRBas  of  Two  Signal  Sources  and  a  Transmitter  of  Opportunity 
(Interferer) 

To  determine  the  CRBas  performance  for  the  selected  conditions,  compare 
Figure  13  (c)  with  Figure  19  (b);  clearly,  a  transmitter  in  the  surveillance  region 
does  not  have  a  big  impact  on  the  CRBas  of  the  moving  signal  source  when  it  is 
in  the  proximity  of  the  stationary  signal  source  in  azimuth  plane.  This  is  true  in 
this  scenario  because  the  transmitter  is  further  apart  from  the  stationary  signal. 
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When  the  CRB  as  of  the  moving  signal  source  is  observed  in  the  azimuth 
region  near  where  the  transmitter  of  opportunity  is  located,  the  CRBas  of  the 
moving  signal  source  is  truly  degraded.  When  the  moving  object  is  at  69°  its 
bound  is  approximately  4  degrees^,  which  is  approximately  equal  to  2°  standard 
variation.  This  value  is  relatively  large,  but  as  the  moving  object  is  within  1°  of 
the  transmitter,  this  increase  in  the  bound  is  crucial. 


Target  SNRs  =  OdB,  Transmitter  SNR  =  50dB  Target  SNRs  =  OdB,  Transmitter  SNR  =  50dB 


(a)  (b) 


(c) 


(d) 
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Target  SNRs  =  OdB,  Transmitter  SNR  =  50dB 


(e) 

Figure  19.  The  CRB  as  of  the  moving  signal  source  with  a  stationary  signal  source 
at  0°  and  a  transmitter  of  opportunity  at  70°  when  SNR  of  the  moving  and 
stationary  signal  sources  are  OdB  and  transmitter  SNR  50dB. 

4.3.2  CRBcm  o£  Two  Signal  Sources  and  a  Transmitter  of  Opportunity 
(Inteiferer) 

The  CRBcm  of  the  moving  object  seems  to  be  a  lower  bound  than  CRBas 
especially  in  the  proximity  of  the  both  stationary  signal  source  and  the 
transmitter  in  azimuth.  As  the  moving  signal  source  moves  away  from  the 
antenna  broadside,  the  bound  increases  due  to  the  array  antenna  performance; 
this  is  shown  in  Section  4.1. 
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(c) 


(d) 


Figure  20.  The  CRBcm  of  the  moving  signal  source  with  a  stationary  signal 


source  at  0°  and  a  transmitter  of  opportunity  at  70°  when  SNR  of  the  moving  and 
stationary  signal  sources  are  OdB,  and  transmitter  SNR  50dB;  and  Start 
Phases/ Doppler  Shifts  of  the  moving,  and  stationary  signal  sources,  transmitter 

are  [2,4,6]/ [8,-2,0]  respectively  . 

Figure  16  (b)  shows  that  the  CRBcm  of  the  moving  signal  source  for  the 
two  signal  case  is  the  same  as  Figure  20  (b).  Figure  20  (b)  shows,  the  CRBcm  of 
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the  moving  object  near  the  stationary  object  when  there  is  also  a  stationary  object 
present  out  at  70°.  This  is  because  of  the  transmitter  being  located  so  far  from  the 
object  of  interest.  In  order  to  show  the  effect  of  a  transmitter  on  the  CRBcm  of  the 
moving  signal  source  in  the  vicinity  of  the  stationary,  the  CRBcm  of  the  moving 
signal  source  is  plotted  in 

Figure  21  considers  the  same  conditions  as  in  Figure  20,  except  that  the 
transmitter  of  opportunity  is  located  at  15°  instead  of  70°. 


Figure  21.  The  CRBcm  of  the  moving  signal  source  with  a  stationary  signal 
source  at  0°  and  a  transmitter  of  opportunity  at  12°  when  SNR  of  the  moving  and 
stationary  signal  sources  are  OdB,  and  transmitter  SNR  50dB;  and  Start 
Phases/ Doppler  Shifts  of  the  moving,  and  stationary  signal  sources,  transmitter 

are  [2,4,6]/ [8,-2,0]  respectively  . 

When  Figure  20  (d)  is  compared  with  Figure  19  (e)  it  is  observed  that  the 
CRBcm  remains  well  below  the  CRBas  when  the  moving  object  is  in  the 
transmitter  region.  Also  comparing  Figure  20  (b)  with  Figure  19  (b)  shows  that 
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bound  of  the  moving  signal  source  is  lower  when  its  estimator  has  the  CM  signal 
constraint.  To  make  a  long  story  short  it  can  be  rephrased  as  follows;  the  CRBcm 
is  a  lower  bound  than  CRB  as  when  the  phase  parameters  of  the  incoming  signals 
have  different  characteristics,  and  the  CRBcm  goes  even  lower  when  the  phase 
parameters  of  each  signal  are  less  correlated. 

4.4  Scenario  4:  Moving  Object  vs.  Different  Transmitter  of 
Opportunity  (Interferer)  Locations 

This  section  demonstrates  the  case  when  there  is  one  transmitter  in  the 
surveillance  region.  The  CRB  as  and  CRBcm  are  presented  below  for  one  moving 
object  signal  whose  DOA  changes  from  -70°  to  70°  for  different  transmitter 
locations.  The  transmitter  location  in  the  azimuth  plane  is  only  changed  from  the 
DOA  spectrum  of  -70°  to  0°  because  the  CRB  is  symmetric  about  the  antenna 
broadside  direction.  This  experiment  should  yield  additional  insight  into  antenna 
placement  and  orientation. 

4.4.1  CRB  AS  of  the  Moving  Object  for  Different  Transmitter  Locations 

As  discussed  earlier,  AS  DOA  estimators  are  more  susceptible  to  increases 
in  the  variance  estimate  when  more  than  one  signal  is  being  received  from  the 
surveillance  region.  This  is  more  pronounced  when  the  signal  of  interest  has 
less  signal  power  compared  with  the  others. 
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Figure  22  show  the  CRBas  of  a  moving  object  when  the  transmitter  is 


located  at  -70°,  -50°,  -30°,  -15°,  0°. 
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Target  SNR  =  OdB.  Transmitter  SNR  =  50dB 


Target  SNR  =  OdB,  Transmitter  SNR  =  50dB 


(i)  0) 

Figure  22.  The  CRB  as  of  the  Moving  Object  with  a  transmitter  of  opportunity  at 
different  locations  [-70°,  -50°,  -30°,  -15°,  0°]  when  SNR  of  the  moving  is  OdB,  and 

transmitter  SNR  50dB. 


The  figures  above  show  that  the  effect  of  the  interfering  transmitter  signal 
on  the  CRBas  decreases  as  the  transmitter  location  shifts  towards  the  antenna 
broadside.  Therefore,  having  the  transmitter  closer  to  the  antenna  broadside 
direction  enables  less  degradation  in  the  bound.  Of  course,  this  results  in  greater 
DOA  estimation  error  in  the  direction  of  the  antenna.  When  the  array  antenna 
orientation  can  be  arbitrarily  chosen,  placing  the  transmitter  out  of  the  coverage 
area  results  in  the  best  accuracy  for  DOA  estimation.  If  it  is  not  possible  to  place 
the  transmitter  out  of  the  coverage  area,  then  the  transmitter  should  be  towards 
an  edge  of  the  array  antenna  where  DOA  estimation  accuracy  is  already 
degraded  due  to  system  configuration.  Although  this  results  in  greater  errors 
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around  the  DOA  of  the  interfering  signal  from  the  transmitter,  the  DOA  for  other 
angular  regions  remains  free  of  interference  and  yields  higher  accuracy 
estimates. 

4.4.2  CRBcm  of  the  Moving  Object  for  Different  Transmitter  Locations 

The  CRBcm  has  been  found  to  be  much  lower  than  the  CRBas  in  the 
presence  of  multiple  signals.  The  results  presented  in  this  section  still  support 
this  conclusion. 

Figure  23  shows  the  CRBcm  of  a  moving  object  when  the  transmitter  is 
located  at  -70°,  -50°,  -30°,  -15°,  0°. 


Start  Phasesl2,4,6)  Doppler[8,-2,0j  Start  PhasesI2.4.6j  Doppleij8,-2.0) 


(a) 


(b) 
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Start  PhasesI2,4.6]  Doppler[8,-2,0] 


True  DOA  (in  degrees) 


(c) 


(d) 


Start  Phases[2,4,6]  Doppler[8 -2,0] 


(e) 


(f) 


77 


CR^Bc^  (degree^)  ^  CRB^^  (degree^) 
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(h) 


Start  Phases[2,4,6]  Doppler[8,*2,OJ 


(i) 


(i) 


Figure  23.  The  CRBcm  of  the  Moving  Object  with  a  transmitter  of  opportunity  at 


different  locations  [-70°,  -50°,  -30°,  -15°,  0°]  when  SNR  of  the  moving  is  OdB,  and 


transmitter  SNR  50dB;  and  Start  Phases/ Doppler  Shifts  of  the  moving,  and 
transmitter  are  [2,6]/ [8,0]  respectively. 


There  is  a  slight  increase  in  CRBcm  when  the 


object  DOA  is  near  the 


transmitter  DOA. 


The  CRBcm  almost  keeps  the  form  of  the  single  object  case  but 
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there  is  a  slight  increase  near  the  transmitter  DOA  when  the  object  is  in  the  same 
azimuth  plane  region. 

In  the  CM  case,  it  is  still  preferable  to  place  the  transmitter  outside  the 
surveillance  region  if  possible.  If  the  transmitter  cannot  be  placed  out  of  the 
surveillance  region,  but  the  antenna  orientation  can  be  varied,  then  the 
transmitter  should  be  placed  towards  the  edge  of  the  array.  By  doing  this,  the 
DOA  region  near  the  antenna  broadside  direction,  which  provides  the  most 
accurate  DOA  estimates,  is  kept  free  from  interference.  It  must  be  noted  that  the 
need  to  change  the  array  antenna  direction  in  the  CM  signal  case  is  not  as  nearly 
pronounced  as  in  the  AS  signal  case. 
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5 .  DISCUSSION 


5.1  Conclusions 

Performance  of  the  Arbitrary  Signal  (AS)  and  Constant  Modulus  (CM) 
signal  DO  A  estimators  were  evaluated  using  the  Cramer-Rao  Bound  (CRB).  The 
CRB  theorem  was  found  to  be  useful  because  it  enables  the  analysis  of  the  error 
variance  for  different  cases.  As  discussed  in  Section  2.3.2,  minimal  error  variance 
is  desirable  and  the  CRB  provides  a  lower  bound  on  this  variance. 

In  the  single  object  case,  both  AS  and  CM  signal  DOA  estimators  achieve 
the  same  bound.  CM  DOA  estimators  use  the  phase  information  of  each 
incoming  signal  to  separate  signal  sources  from  each  other  and  DOA  estimation 
is  done  separately  for  each  signal.  However,  when  only  one  signal  is  present  in 
the  environment,  the  CRBcm  is  identical  to  the  CRB  as. 

In  the  case  of  multiple  signals,  the  CRBcm  is  a  lower  bound  around  the 
region  where  the  angular  separation  between  two  signals  in  azimuth  gets 
smaller,  even  when  one  of  the  signals  is  strong  interference  from  a  transmitter. 
The  CRBcm  is  found  to  be  slightly  above  the  CRB  for  the  single  object  case, 
including  the  case  with  the  interfering  transmitter.  On  the  other  hand,  the  CRB  as 
for  Arbitrary  Signals  increases  greatly  when  the  angular  separation  between  two 
signals  becomes  small,  and  CRBas  is  also  affected  dramatically  by  interfering 
signals  with  large  signal  power,  (such  as  the  interfering  transmitter  signal). 
Hence,  this  research  demonstrated  that,  in  general,  the  bounds  of  the  CM  DOA 
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estimators  are  lower  than  the  bounds  for  AS  DOA  estimators;  CM  DOA 


estimators  are  declared  to  have  better  performance. 

As  mentioned  in  Section  2.3,  the  CRB  shows  the  performance  of  a 
Minimum  Variance  Unbiased  (MVUB)  estimator.  In  spectral  estimation,  the 
search  for  a  MVUB  estimator  is  a  challenging  task.  Since  the  bounds  for  CM 
DOA  estimators  show  good  characteristics  in  scenarios  closely  representing  real- 
world  applications,  the  CM  DOA  estimators  prove  to  be  promising.  Hence,  in 
the  search  for  the  MVUB,  giving  consideration  to  CM  Signal  DOA  estimators  is 
highly  encouraged. 

5.2  Future  Work 

In  this  thesis  a  single  Uniform  Linear  Array  (ULA)  configuration  was 
studied  and  there  was  only  one  DOA  estimate  obtained  from  the  array.  Dividing 
the  array  into  smaller  sub-arrays  results  in  degraded  estimation  performance  [6]. 
On  the  other  hand,  combining  all  sub-array  estimates  and  creating  mutual 
coupling  between  the  sub-arrays  resulted  in  a  new  estimator  with  different 
performance.  The  change  in  overall  DOA  estimation  performance  should  be 
studied  and  its  performance  compared  to  a  non-divided  array.  As  demonstrated 
in  this  thesis,  the  CRB  provides  an  excellent  tool  for  such  a  study. 

As  was  demonstrated  in  Sections  4.2  and  4.3,  the  DOA  CRB  of  a  signal 
when  other  signals  exist  in  the  environment  show  different  characteristics.  The 
CRB  was  found  to  increase  when  the  DOA  of  a  signal  gets  close  to  the  DOA  of 
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other  signals.  This  is  especially  true  when  there  is  correlation  between  the 
signals;  the  increase  in  the  CRB  becomes  more  significant.  Correlation  effects 
between  two  signals  was  implemented  in  this  research  using  a  correlation 
coefficient  and  the  effect  of  the  correlation  level  between  the  Arbitrary  Signals  on 
their  CRB  was  shown.  However,  modeling  correlation  between  CM  signals  is 
not  as  simple.  In  the  thesis  it  was  mentioned  that  worst  case  for  CM  signals 
occurs  when  the  phase  parameters  of  two  signals  coming  from  the  same 
direction  are  exactly  the  same.  This  is  the  case  when  the  phases  of  two  signals 
are  fully  correlated.  The  changes  in  the  CRBcm  with  respect  to  correlation  level 
of  the  two  signals  with  close  DOA  angles  should  be  studied  by  parameterizing 
the  correlation  level  of  the  CM  signals.  This  would  possibly  require  a  new 
derivation  of  the  CRB  for  correlated  CM  signals. 

One  assumption  that  simplified  the  extent  of  the  research  is  that  the 
objects  are  zero  altitude  objects  with  DOA  only  changing  in  azimuth.  Adding 
elevation  to  the  individual  objects  and  including  their  altitude  in  the  calculations 
will  reflect  results  that  are  more  suited  to  real-world  applications.  Bounds  for 
DOA  error  in  elevation  and  azimuth  can  be  achieved  by  using  a  planer  array 
antenna  model. 

The  numerical  calculations  of  the  CRB  are  based  on  the  system  described 
in  Section  1.5,  which  is  composed  of  a  single  ULA  antenna  configuration.  As 
mentioned  above,  this  can  easily  be  extended  to  sets  of  smaller  ULA  antennas, 
although  the  composite  estimate  becomes  somewhat  more  complicated. 
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However,  ULA  antenna  configurations  are  only  one  of  many.  Cramer-Rao 
Bounds  for  other  types  of  array  configurations  can  also  be  determined  and  the 
advantages/ disadvantages  of  different  array  types,  shapes,  radiating  elements 
and  element  spacing  can  be  studied. 
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Appendix  A 


PRF  1 :  E  e^(t)e(t)e^(s)e(s)  = 


fortes 
m(m  +  l)a^  fort  =  s 


When  t  /  s,  e(t)  and  e(s)  are  statistically  independent  from  each  other 
hence  the  expected  value  of  their  product  can  be  written  as  the  products  of  their 
expected  values. 

Et"(t)e(t)e"(s)e(s)]=E|e"(t)e(t)]E|e"(s)e(s)]=yE[e,’(t)e,(t)]x  jE[e|’(s)e,(s)] 


Z-xE 

i=l  1=1 


2  2 

a  =  m  a 


When  t  =  s. 


E^"(t)e(t)e“(t)e(t)]=  E[(e"(t)e(t)t]  =  E[(g^(t)i(t)+ " 

=  E[fe^(t)e(t))n  +  2E[(i^(t)i(t))(e^(t)e(t))]+Er(e^(t)e 


Note  that  the  real  and  imaginary  parts  are  statistically  independent. 

E[e“  (t) e(t) e“  (t) e(t)]  =  (t)  e(t)f  +  2E[e^  (t)  e(t)]E[(e ^  (t)  e (t))]+  E  (e ^  (t)  e ' 


E  e^(t)e(t)  =  E  e^(t)e( 


E^W^t);  =E(e^(t)ei 


E[e«(t)e(t)e«(t)e(t)]=2Erfe^(t) 


)e(ttfh^ 


Erfe^(t)e(t)n  =  Ef|:e.^(t)l  =Erje.^(t)x|;e.^(t)l  =  X|;E[e.^(t)e,^(t)]+XEk^^ 


J  (t)]+  m  E[e/  (t) 
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E  (e^(t)e(t))  =  (m -l)m^  + =  (m  + 2)m^ 

^  E  e”  (O^W (O^W]  =  '''  ^2  ~  '*' 

PRF  2:  E  e^  (t)  e(t)  e”  (s)  =0 

When  t  #  s,  they  are  independent,  and  because  the  noise  has  zero  mean, 
the  proof  is  strait  forward. 

E  e”  (t)  e(t)  e”  (s)  =E  e”(t)  e(t)  E  (s)  =0 

When  t  =  s,  the  expectation  becomes  the  third  order  moment  of  a  Gaussian 
random  variable,  which  is  zero. 


This  equality  can  also  be  written  as 


HG+HG=I 

HG+HG=0 


A 

which  must  hold  true  since  G  =  H“',  and 

I  =  HG  =  (h  +  jH)(G  +  jG)=  (hG  -  Hg)+  j(HG  +  Hg)  . 


85 


Appendix  B 


CRB^s(0r)  =  -  ^Re  X“(t)D^  I-A(a“a)  ‘a“  DX(t)  is  given  in  (52) 

2  t=i  L  L  J  A  ) 


Then  the  Jij  can  be  shown  as. 


Re  [l  -A(a''a)‘a“1  d(0Rj)  ^x“(t)xj(t) 

L  >-  -I  t=i 


A  r  1  N 

P  =  E  x(t)x“(t)J  and  by  definition  Substituting 

N^oo  N  t=l 

Pij  into  Jij  and  rewriting  Fisher  Information  matrix  gives. 


-Rejo”  I-A(a“a)'‘a“  D.*NP^  . 


Then  the  CRB  as  is  obtained  as  given  in  Equation  (53). 


CRB^,  (0  J  =  ^  (Re[D»  [l  -  a(a“  a)-‘ A»  ]d.*  P"  ] )‘‘ 
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Appendix  C 

P  =  e[xx“] 

When  n-  2, 

E[xf]  EfxjxJ 
LE[xiX2]  e[x^] 

For  convenience  let's  assume  thatEj^x^  j=l,  and  Ej^Xj  j  =  a  e[x^ J  so  P  can 
be  written  as, 

p^r  1  EfxjXj] 

LE[xiX2]  a 

E  f  X  X  ]  j — 

Correlation  coefficient,  p  is  equal  to  Pj2  = - ,  Since  =  vaa^  ,  and 

a  a  2  1 

x,  X2 

=1;  e[xjX2]=  Pi2a^  =Pj2Va.  Substituting  this  into  P  gives. 


1  P[2  Vtt 

Pj2  Vtt  a 
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Appendix  D 

This  appendix  shows  the  CRB  for  a  single  signal  over  a  DOA  range  of  70° 
to  -70°,  along  with  effects  of  different  SNR  levels.  In  Section  4.1,  the  CRB  for  a 
single  signal  is  shown  when  SNR  is  OdB.  Here  the  CRB  of  a  single  signal  is 
plotted  for  the  SNR  levels  of  10,  20,  -lOdBs.  In  numerical  representation  of  the 
SNR  this  corresponds  to  the  increase  of  SNR  from  OdB  by  10,  100  and  1/10 
respectively.  There  is  an  inverse  relation  between  the  SNR  and  CRB.  When  the 
SNR  is  increased  by  factors  of  10,  100,  and  1/10;  CRB  is  reduced  by  factors  of 
1/10,  1/100,  10.  This  can  be  observed  in  the  plots  below.  Please  note  that  the 
shape  remains  the  same  but  the  magnitude  on  the  axis  where  CRB  is  shown 
changes  by  the  factors  that  have  been  mentioned  above. 


SNR  =  -10dB 
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SNR  =  10dB 


x10"' 


SNR  =  20dB 
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Appendix  E 


CRBas  of  the  moving  and  the  stationary  signal  sources  are  given  in  the 
figures  below.  As  the  correlation  between  the  signal  sources  increases,  the 
CRBas  increase  dramatically  when  two  signal  sources  are  close  to  one  another  in 
the  azimuth  plane.  When  the  signal  sources  are  further  apart  from  each  other  in 
azimuth,  then  CRBas  is  still  higher  for  the  highly  correlated  signals  case  but  it  is 
not  as  pronounced  as  when  angular  signal  separation  of  the  signals  is 
comparatively  less. 

Moving  Target  SNR  =  OdB 


0.012 


0.5 

Correlation  °  ^ 


20  35  50 

True  DOA  (in  degrees) 
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CRB^g  (degree^) 


Moving  Target  SNR  =  OdB 


True  DOA  (in  degrees) 


The  CRB  AS  of  the  moving  signal  source  when  there  is  a  stationary  signal  source 
located  at  antenna  broadside  with  OdB  SNR  values  for  both. 
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stationary  Target  SNR  =  OdB 


0.01 


c<^.008 

0) 

O) 

^0.006 

(n 

< 

CQ 

g  0.004 
0.002 


Correlation 


Angular  Separation  (in  degrees) 


The  CRB  AS  of  the  stationary  signal  source  located  at  the  antenna  broadside  when 


there  is  a  moving  signal  source  getting  closer  to  the  stationary  signal  source  in 


azimuth  plane  with  OdB  SNR  values  for  both. 


92 


BIBLIOGRAPHY 


[1]  Skolnik,  Merrill  I.  Introduction  to  Radar  Systems  (3rd  Edition).  New  York: 
McGraw  Hill  Co.,  2001. 

[2]  Howland,  Paul  E.  “A  Review  of  Passive  Sensor  Technology,”  Surveillance  Branch, 
NATO  C3  Agency,  Netherlands,  14  Sep  2000. 

[3]  Willis,  Nicholas  J.  Bistatic  Radar  (2nd  Edition).  Maryland:  Technology  Service 
Corporation,  1995. 

[4]  Kay,  Steven  M.  Eundamentals  of  Statistical  Signal  Processing.  NJ:  Prentice  Hall, 
1993. 

[5]  Scharf,  Eouis  E.  Statistical  Signal  Processing  Detection,  Estimation,  and  Time 
Series  Analysis.  Massachusetts:  Addison-Wesley  Publishing  Co.  Inc.,  1991. 

[6]  Stoica,  Petre,  and  Arye  Nehorai.  ‘MUSIC,  Maximum  Eikelihood,  and  Cramer-Rao 
Bound,”  IEEE  Transactions  on  Acoustics,  Speech,  and  Signal  Processing,  37:  720- 
741  (May  1989) 

[7]  Van  der  Veen,  Alle-Jan,  and  Arogyaswami  Paulraj.  “An  Analytical  Constant 
Modulus  Algorithm,”  IEEE  Transactions  on  Signal  Processing,  44:  (May  1996). 

[8]  Eaves,  J.  E.  and  E.  K  Reedy.  Principles  of  Modern  Radar.  New  York:  Van 
Nostrand  Reinhold,  1987. 

[9]  Maas,  S.  A.  “Microwave  Mixers  in  the  90s,”  Microwave  Journal,  61-72  Mar  1990. 

[10]  Guerlac,  H.  E.  Radar  in  World  War  II,  Vols.  I  and  II.  New  York: 
Tomask/ American  Institute  of  Physics,  1987. 

[11]  Skolnik,  M.  I.  “Hn  analysis  of  Bistatic  Radar,""  IRE  Transactions  on  Aerospace  and 
Navigational  Elektroniks,  19-27,  Mar  1961. 

[  12]  Price,  A.  Instruments  of  Darkness  -  The  Struggle  For  Radar  Suppremacy.  William 
Kimber  and  Co.  Etd.,  1967. 

[13]  Glaser,  J.  I.  ''Some  Results  in  the  Bistatic  Cross  Section  of  Complex  Objects,” 
Proceedings  of  the  IEEE,  11,  5:  639-648  (May  1989). 

[14]  Willis,  Nicholas  J.  “PCE  Systems.”  Address  to  Air  Eorce  Institute  of  Technology 
(AEIT)  Students  and  Eaculty  Members.  AEIT,  Wright-Patterson  APB,  Ohio.  June 
2001. 


[15]  Veen,  Alle  J.  “Blind  Souree  Separation  Based  On  Combined  Direetion  Finding  and 
Constant  Modulus  Properties,”  Delft  University  of  Teehnology,  Netherlands. 

[16]  Howland,  Paul  E.  Television  Based  Bistatic  Radar,  PhD  thesis.  Sehool  of 
Eleetronie  and  Eleetrieal  Engineering,  University  of  Birmingham,  England,  Sep 
1997  (B152TT). 

[17]  Stimson,  George  W.  Introduetion  to  Airborne  Radar  Edition).  New  Jersey: 
Seiteeh  Publishing  Ine.,  1998. 

[18]  Balanis,  Constantine  A.  Antenna  Theory  Analysis  and  Design  (2"‘*  Edition).  New 
YorkiJohn  Wiley  &  Sons  Ine.,  1997. 

[19]  Skolnik,  Merrill  1.  Radar  Handbook  (2"‘*  Edition).  New  York:  MeGraw  Hill  Co., 
1990. 

[20]  Eraneos,  Joseph  M.,  and  Benjamin  Eriedlander.  “Bounds  for  Estimation  of 
Complex  Exponentials  in  Unknown  Colored  Noise,”  IEEE  Transactions  on  Signal 
Processing,  43:  2176-  2185  (September  1995). 

[21]  Stark,  Henry,  and  John  W.  Woods.  Probability,  Random  Proeesses,  and  Estimation 
Theory  for  Engineers  (2"^^  Edition).  New  Jersey:  Prentiee  Hall  Ine.,  1994. 

[22]  Gooeh  R.,  and  J.  Eundell.  “The  CM  Array:  An  Adaptive  Beamformer  for  Constant 
Modulus  Signals,”  IEEE  International  Conference  on  Acoustics,  Speech,  and 
Signal  Processing,  2523-2526  (May  1986). 

[23]  Treieher,  J.R.  and  M.G.  Earimore.  “New  proeessing  teehniques  based  on  eonstant 
modulus  adaptive  algoritm,”  IEEE  Transactions  on  Acoustics,  Speech,  Signal 
Processing,  33:  420-431  (April  1985). 

[24]  Leshem,  Amir  and  Alle  J.  Van  der  Veen.  “Direetion-of- Arrival  Estimation  for 
Constant  Modulus  Signals,”  IEEE  Transactions  on  Signal  Processing,  47:3125- 
3130  (November  1999). 

[25]  Sadler,  Brian  M.,  and  Riehard  J.  Koziek.  “Bounds  on  Bearing  and  Symbol 
Estimation  with  Side  Information,”  IEEE  Transactions  on  Signal  Processing,  49: 
822-834  (April  2001). 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  074-0188 

The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  the  collection  of  information,  including 
suggestions  for  reducing  this  burden  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway, 

Suite  1 204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  an  penalty  for  failing  to  comply  with  a  collection  of 
information  if  it  does  not  display  a  currently  valid  0MB  control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 

1.  REPORT  DATE  (DD-MM-YYYY) 
15-03-2002 

2.  REPORT  TYPE 

Master’s  Thesis 

3.  DATES  COVERED  (From  -  To) 

Jun  2001  -  Mar  2002 

4.  TITLE  AND  SUBTITLE 

STATISTICAL  ERROR  ANALYSIS  OF  A  DOA  ESTIMATOR  FOR  A  PCL 
SYSTEM  USING  THE  CRAMER-RAO  BOUND  THEOREM 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 


6.  AUTHOR(S) 

Say,  Kerim,  U‘  Lt.,  TUAF 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAMES(S)  AND  ADDRESS(S) 

Air  Force  Institute  of  Technology 

Graduate  School  of  Engineering  and  Management  (AFIT/EN) 
2950  P  Street,  Building  640 
WPAFB  OH  45433-7765 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

NATO  C3  Agency 
Attn:  Dr.  Paul  Howland 

P.OBoxl74  Telephone  No:  -^3 1-70  314-2476 

The  Netherlands  e-mail:  Paul.Howland@nc3a.nato.int 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED. 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

AFIT/GE/ENG/02M-23 


10.  SPONSOR/MONITOR’S  ACRONYM(S) 


11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 


14.  ABSTRACT 

Direction  of  Arrival  (DOA)  estimation  of  signals  has  been  a  popular  research  area  in  Signal  Processing.  DOA  estimation  also  has 
a  significant  role  in  the  object  location  process  of  Passive  Coherent  Location  (PCL)  systems.  PCL  systems  have  been  in  open 
literature  since  1986  and  their  applications  are  not  as  clearly  understood  as  the  DOA  estimation  problem.  However,  they  are  the 
focus  of  many  current  research  efforts  and  show  much  promise. 

The  purpose  of  this  research  is  to  analyze  the  DOA  estimation  errors  in  a  PCL  system.  The  performance  of  DOA  estimators  is 
studied  using  the  Cramer-Rao  Bound  (CRB)  Theorem.  The  CRB  provides  a  lower  bound  on  the  variance  of  unbiased  DOA 
estimators.  Since  variance  is  a  desirable  property  for  measuring  the  accuracy  of  an  estimator,  the  CRB  gives  a  good  indication  about 
the  performance  of  an  estimator.  Previous  DOA  estimators  configured  with  array  antennas  used  the  array  antenna  manifold,  or  the 
properties  of  the  array  antenna  structure,  to  estimate  signal  DOA.  Conventional  DOA  estimators  use  arbitrary  signal  (AS)  structures. 
Constant  Modulus  (CM)  DOA  estimators  restrict  the  input  signals  to  a  family  of  constant  envelope  signals,  and  when  there  are 
multiple  signals  in  the  environment,  CM  DOA  estimators  are  able  to  separate  signals  from  each  other  using  the  CM  signal  property. 
CM  estimators  then  estimate  the  DOA  for  each  signal  individually. 

This  research  compares  the  CRB  for  AS  and  CM  DOA  estimators  for  a  selected  system.  The  CRB  is  also  computed  for  this 
system  when  single  and  multiple  and  moving  objects  are  present.  The  CRB  as  and  CRB^m  are  found  to  be  different  for  the  multiple 
signal  case  and  moving  object  cases. _ 


15.  SUBJECT  TERMS 

Passive  Radar,  Passive  Coherent  Location  System,  Passive  Sensor  Location  System,  Bistatic  Radar,  Multistatic  Radar,  Cramer- 
Rao  Bound,  Direction  of  Arrival,  DOA  Estimator,  Array  Antennas,  Statistical  Error  Analysis. 


X 


standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  Z39-18 


